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1.  Introduction 

The  intent  of  this  article  is  to  study  the  behavior  of  solutions  (v(x,t),u(x,t))  of  the 


system  of  differential  equations 

Vt  =  (T  x 

(U) 

ut  =  vx 

(1.2) 

where 

a  =  r(ti)t>£ 

(1.3) 

with  t(u)  a  smooth  function  satisfying 

r(u)  >  0  ,  t'(u )  <  0 

(1.4) 

and  n  a  positive  parameter.  Equations  (1.1  -  1.3)  give  rise  to  a  coupled  system  of  par¬ 
tial  differential  equations  in  one  space  dimension.  They  are  supplemented  with  initial 
conditions 

u(z,0)  =  uo(z)  ,  u(x,  0)  =  u0(x) ,  (1.5) 

and,  as  a  consequence, 

<t(i,0)  =  a0(x)  :=  t(u0(x))v£x(x)  ,  (1.6) 

and  with  boundary  conditions  that  are  discussed  later. 

To  gain  some  perspective  on  the  problem,  note  that,  if  n  =  0,  (1.1  -  1.3)  leads  to  the 
pair  of  conservation  laws 

v t  =  t(u)z 

V  (1.7) 

ut=vx. 

If  t'(u )  >  0  then  (1.7)  is  hyperbolic;  however,  under  (1.4),  the  system  (1.7)  is  elliptic  and 
the  initial  value  problem  is  ill-posed.  Nevertheless,  it  admits  an  interesting  class  of  special 
solutions 

v(x,  t)  =  x 
u(x,  t)  =  t  +  Uq  , 


(l.S) 


where  uo  is  an  arbitrary  constant.  Equations  (1.1  -  1.3)  with  n  >  0  can  be  thought  as  a  ~ 
particular  regularization  of  (1.7).  ~ 
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A  motivation  for  studying  this  problem  stems  from  a  program  of  understanding  the 
phenomenon  of  shear  band  formation  at  high  strain  rates.  Shear  bands  are  narrow  regions 
of  intensely  concentrated  shearing  deformation  that  are  observed  during  the  plastic  de¬ 
formation  of  many  materials.  The  occurrence  of  shear  bands  is  typically  associated  with 
strain  softening  type  response,  past  a  critical  strain,  of  the  measured  average  shear  stress 
£(t)  versus  the  measured  average  shear  strain#Z7(<);  that  is,  £  =  r(Lr),  where  r(-)  is  in¬ 
creasing  up  to  a  certain  critical  strain  and  decreasing  thereafter.  Various  mechanisms  and 
associated  continuum  thermomechanics  models,  often  depending  on  the  particular  context, 
have  been  proposed  for  the  explanation  of  shear  bands^(see  Shawki  and  Clifton  [15]  for  an 
excellent  survey  of  the  related  literature),*An  underlying  common  feature  of  several  models 
is  that  they  are  regularizations  of  an  ill-posed  problem,  or  that  the  associated  linearized 
problem  exhibits  growth  of  high  frequency  modes,  *'  •  ■  1  ■■  1  v 

The  model  employed  here  describes  the  plastic  shearing  of  an  infinite  plate  of  unit 
thickness  subjected  to  either  prescribed  tractions  or  prescribed  velocities  at  the  boundaries. 
In  this  framework,  u(r,  t )  describes  the  velocity  field  in  the  shearing  direction,  <r(x,  t)  stands 
for  the  shear  stress  and  u(x,t)  for  the  plastic  shear  strain.  Equation  (1.1)  describes  the 
balance  of  linear  momentum,  while  (1.2)  is  a  kinematic  compatibility  relation  (note  that 
elastic  effects  are  neglected);  (1.1)  and  (1.2)  are  taken  over  (r,t)e[0,l]  x  {f  >  0}  and  are 
supplemented  with  the  boundary  conditions 

-  f- 

in  case  the  shearing  deformation  is  caused  by  prescribed  tractions  at  the  boundaries,  or 

v(0,t)  =  0,  u(l,<)  =  1  ,  f>0,  (l-9)r 

in  the  case  of  prescribed  velocities.  The  constitutive  law  (1.3)  is  appropriate  for  a  material 
exhibiting  strain  softening,  as  manifested  in  (1.4),  and  strain  rate  sensitivity,  the  strength 
of  which  is  measured  by  the  parameter  n.  Our  objective  is  to  use  (1.1  -  1.6),  (1.9)  as  a 
test  problem  to  analyze  the  competition  between  the  destabilizing  effect  of  strain  softening 
versus  the  stabilizing  effect  of  strain  rate  dependence. 

Technically,  the  model  (1.1  - 1.3)  belongs  to  the  class  of  isothermal  viscoelasticity  of  the 
rate  type  (for  general  information  on  the  mathematical  theory  of  viscoelasticity  the  reader 
is  referred  to  Renardy,  Hrusa  and  Nohel  [14]).  Metals,  in  general,  exhibit  strain  hardening 
in  isothermal  deformations.  However,  an  increase  in  temperature  causes  a  decrease  in  the 
yield  stress,  so  that  in  an  adiabatic  deformation  the  combined  effect  of  strain  hardening 
and  thermal  softening  may  deliver  eventually  a  net  softening.  Thus,  although  (1.1  -  1.3)  is 
a  model  in  the  framework  of  isothermal  mechanical  theories,  thermal  effects  are  implicitly 
taken  into  account  through  the  hypothesis  of  strain  softening  .  One  of  our  goals  is  to 
reveal  similarities  in  the  structure  and  predictions  of  (1.1  -  1.3)  as  compared  to  related 
models  incorporating  thermal  effects  that  have  been  studied  recently  in  the  mathematical 
literature  [7,  17,  19,  2,  1], 

We  emphasize  that  the  spirit  of  this  study  is  not  to  recover  solutions  of  ( 1.7)-(  1.4)  as 
n  — »  0  limits  of  solutions  of  (1.1  -  1.4).  Rather  the  rational  here  is  the  converse.  Because 
of  the  inherent  instability  induced  by  strain  softening  ,  it  has  been  postulated  that  higher 


cr(0,t)  =  1 ,  <r(l,t)  =  1  ,  <>0, 


(l-9)< 


order  effects,  such  as  strain-rate  dependence,  play  an  important  role  and  cannot  be  ignored 
(cf.  [11,  21]).  Apart  from  some  previous  investigations  using  (1.3)  for  n  =  1  [18,  2],  other 
types  of  rate  dependent  constitutive  relations  have  been  used  to  analyze  shear  bands  (e.g. 
Wu  and  Freund  [21]),  as  well  as  strain-gradient  dependent  constitutive  laws  (e.g.  Coleman 
and  Hodgdon  [5]).  There  is  a  very  extensive  mechanics  literature  on  the  subject  and  the 
reader  is  referred  to  [21,  11,  20,  15]  and  references  therein. 

From  an  analysis  point  of  view,  equations  (1.1  -  1.3)  give  rise  to  a  coupled  system 
consisting  of  a  parabolic  equation  in  v  coupled  through  the  diffusion  coeffient  with  (1.2)  (cf. 
(3.1)).  As  the  material  is  being  sheared,  under  the  effect  of  (1.9),  the  diffusion  coefficient 
is  decreasing.  It  is  conceivable,  that  if  the  decrease  is  too  rapid  and/or  nonuniform  in 
the  space  variable,  the  diffusion  may  not  be  able  to  stabilize  the  process.  To  analyze 
this  competition,  it  is  helpful  to  recast  (1.1  -  1.3)  into  an  equivalent  formulation  of  a 
reaction-diffusion  system  (cf.  (3.4  -  3.5)). 

In  Section  2  we  pursue  an  existence  theory  of  classical  solutions  for  a  coupled  system 
of  partial  differential  equations  (cf.  (2.1  -  2.3))  that  includes  (1.1  -  1.3).  This  system  also 
includes  certain  more  general  models  in  viscoelasticity  with  internal  variables,  as  well  as 
some  models  incorporating  thermal  effects  that  are  used  for  the  analysis  of  shear  bands 
[15,  19].  Motivated  by  the  problems  under  consideration,  the  main  objective  is  to  identify 
a  minimal  set  of  a-priori  estimates  sufficient  for  continuation  of  solutions.  The  existence 
theory  is  done  in  Schauder  spaces  and  the  main  ingredient  is  an  application  of  the  Lerav- 
Schauder  fixed  point  theorem.  The  results  are  summarized  in  Theorems  2.4  and  2.5.  For 
existence  theories  of  weak  solutions  in  structurally  related  systems  the  reader  is  referred 
to  Charalambakis  and  Murat  [3]  and  Nohel  et  al.  [12]. 

In  Section  3  we  take  up  the  problem  (V)s  consisting  of  (1.1  -  l.G)  with  stress  boundary 
conditions  (1.9)5.  Using  the  results  of  Section  2  together  with  the  special  structure  of  the 
system,  it  is  shown  in  Theorem  3.2  that  solutions  of  ( V)s  are  globally  defined  if  and  only  if 
the  integral  r(£)"  d£  diverges.  Moreover,  the  evolution  of  solutions  of  ('P)s  is  studied 
under  various  assumptions  for  the  constitutive  function  r(u).  Below,  we  summarize  the 
outcome  of  the  analysis  for  the  special  case  of  a  power  law 


with  parameters  m,  n  positive.  The  parameter  region  is  decomposed  into  three  distinct 
subregions  0<-^<-~,|<-^<l  and  —  >  1,  across  which  the  response  changes 
drastically: 

(i)  In  the  region  0  <  ~  <  4  solutions  of  (V)s  are  globally  defined  and,  as  t  — ►  oc,  the  shear' 
stress  a(x,t)  is  attracted  to  the  constant  state  (7  =  1  while  u(.r,t)  behaves  asymptotically 
as  a  function  of  time. 

(ii)  In  the  region  4  <  ^  <  1  the  constant  state  <7  =  1  loses  its  stability  and  nonuniformities 
in  the  strain  persist  for  all  times. 

(iii)  Finally,  in  the  region  ^  >  1,  u(x,t)  becomes  infinite  in  finite  time. 

For  ^  >  2'  we  exhibit  initial  data  for  which  the  corresponding  u(x,t)  develops  nonunifor¬ 
mities  around  x  =  0  and  x  =  1  and  looks  like  two  shear  bands  located  at  the  boundaries. 
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The  analysis  of  Section  3  is  effected  by  means  of  comparison  principles  for  (3.4)  and  ener°w 
estimates  for  (V)s.  6 

For  a  power  law  (1.10)  the  system  (1.1  -  1.2)  is  invariant  under  a  family  of  scaling 
transformations.  In  Section  4  we  take  up  the  problem  (V)y  consisting  of  (1.1  -  1.2),  (1.10) 
and  velocity  boundary  conditions  (1.9)  v,  and  introduce  a  change  of  variables  motivated 
by  the  scaling  property.  The  resulting  system  (4.20  -  4.22)  admits  positively  invariant 
rectangle^  of  arbitrary  size.  Using  this  observation  together  with  energy  estimates  for 
{V)v,  it  is  shown  m  Theorem  4.1  that,  if  m  <  min{n ,  1},  every  solution  of  (V)v  converges 
to  the  uniform  shearing  solution  (1.8),  as  t  -+  oo. 


2.  Existence  Theory  and  Regularizing  Effect  for  a  Coupled  System 

We  consider  the  initial-boundary  value  problem  consisting  of  the  system  of  quasilinear 
partial  differential  equations 

dtv  =  dzw  (2.1) 

dtu  =  f(x,u,vx)  (2.2) 

for  ( x,t)eQT  :  =  (0,1)  x  (0,T],  T  >  0,  where 


with  boundary  conditions 


w  =  <f(x,u,vz) , 


(2.3) 


v(0,t)  =  v(l,t)  =  0,  0  <t<T,  (2.4)v- 

or 

w{0,t)  =  w(l,t)  =  0,  0<t<T,  (2.4)5 

and  initial  conditions 


u(i,0)  =  u0(j),  u(j,0)  =  u0(x),  0  <  x  <  1 ;  (2.5) 

as  a  consequence  of  (2.3)  and  (2.5), 

w(x,  0)  =  u>o(x)  :=  <p(x,  u0(x),  v0z(x )) .  (2.6) 

The  functions  v(x,  t),  w(x,t)  are  real  valued,  while  u(x,  t)  stands  for  an  Revalued  func¬ 
tion,  all  defined  on  Qt  —  [0, 1]  x  [0,  T\.  The  given  functions  /(. r,  p,  q)  :  [0,1]  • .  Its  xR  — ►  RN 
and  <p(x,p,  q)  :  [0, 1]  x  RN  x  R  -+  R  are  assumed  to  be  smooth  with  respect,  to  all  their 
arguments  (the  hypothesis  /  and  of  class  C 2  suffices  for  all  that  follows).  In  addition 
for  each  fixed  (x,p)  the  function  <p(x,p,  •)  is  assumed  to  be  strictly  increasing  and  thus 

invertible.  Let  i/’(x,p,r)  :  [0,1]  x  RN  x  R  — ►  R  be  the  inverse  function.  Inverting  (2.3) 
yields 

vz  =  V’(.t,u,uj)  .  (2j) 
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We  seek  solutions  (v(x,t),u(x,t))  of  (2.1  -  2.5)  defined  on  Qt,  T  >  0.  Our  specific 
goals  are  to  identify  a  minimal  set  of  a-priori  estimates  that  guarantee  existence  and 
continuation  of  solutions  up  to  time  T  >  0,  and  to  study  the  regularizing  effect  that  the 
parabolic  equation  (2.1),  (2.3)  exerts  on  solutions. 

To  this  end  it  is  expedient  to  state  an  alternative  formulation  of  the  problem.  The 
initial-boundary  value  problem  (2.1  -  2.5)  is  formally  equivalent  to  the  system  of  reaction- 
diffusion  equations 

dtw  —  a(x,u,w)d\w  +  6(x,  u,  w)  (2.8) 

dtu  =  g(x,u,w )  (2.9) 

with  boundary  conditions 

ti>*(0,<)  =  «>*(1,<)  =  0,  0  <  f  <  T ,  (2.10)v 

or 

w(0,t)  =  w(l,t)  =  0,  0  <  t  <  T ,  (2.10)5 

and  initial  conditions 


w(x,  0)  =  w0(x) ,  u(x,  0)  =  u0(x) ,  0  <  x  <  1 ;  (2.11) 


where 

a(x,  u,w)  —  tfq(x,U,  tp(x,  u.  w )) 

6(x,  u,  w)  =  {ipp  ■  /)(x,  u,  t/>(x,  u,  w ))  (2.12) 

g(x,u,  w)  =  /(x,u,0(x,it,iw)). 


Indeed,  given  any  sufficiently  smooth  solution  (v(x,  t),  u(x,  t))  of  (2.1  -  2.5)  the  pair 
(w(x,  t),  u(x,  t))  satisfies  (2.8  -  2.11)  as  follows:  Differentiating  (2.3)  with  respect  to  t  and 
using  (2.1),  (2.2),  (2.7)  and  (2.12)  leads  to  (2.8);  (2.2)  and  (2.7)  yield  (2.9);  the  rest  are 
clear.  Conversely,  if  (w(x,t),u(x,t))  is  a  classical  solution  of  (2.8  -  2.12),  define  a  function 
v(x,t)  on  Qt  such  that 


vz  =  t J>(x,u,w) 
vt  ^  wx 


(2.13) 


and  v(x,0)  =  vq(x),  0  <  x  <  1.  Since  ip(x,p,  ■)  is  the  inverse  function  of  <p(x,p,  •),  the  list 
of  relations 

<p{x,p,ij>(x,p,r))  =  r 

V??(x,p,^(x,p,r))0r(x,p,r)  =  1  (2.14) 

<p?(x,p,t/>(x,p,r))  +  <p9(x,p.  t/>(x,p,  r))ipp(x,p,  r)  =  0 

holds,  and  the  compatibility  of  (2.13)  amounts  to  (2.8)  via  (2.9)  and  (2.12).  Moreover. 
(u(x,  t),  u(x,  t))  satisfies  (2.1  -  2.5). 

Our  strategy  is  to  first  prove  an  existence  theorem  for  classical  solutions  of  (2.S  - 
2.11),  in  Schauder  spaces,  using  the  Leray-Schauder  fixed  point  theorem  ([10]).  This,  in 
turn,  yields  an  existence  theorem  for  the  equivalent  system  (2.1  -  2.5)  provided  the  initial 
data  are  sufficiently  smooth.  The  smoothness  assumptions  are  then  relaxed  by  means  of 
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density  arguments.  In  the  sequel  j-|  will  stand  for  both  the  absolute  value  and  the  Euclidean 
norm  in  RN .  Also,  ||  •  {{0,13/2  (I  •  U)  will  denote  the  usual  Schauder  norms  (cf.  [8,9])  in 

/2(Qt)  (Ca[0, 1])  or  \CP'PJ2(Qt)\  N  ([C^O,  1]] N)-  The  meaning  of  these  symbols  will 
be  apparent  from  the  context. 

The  possibility  that  (2.1  -  2.5)  (or  (2.8  -2.11))  admit  globally  defined  solutions  for 
general  nonlinear  functions  can  be  ruled  out  by  considering  special  cases  when  the  system 
decouples.  To  ensure  global  solvability  one  could  place  certain  growth  restrictions  on  the 
functions  /  and  9  (or  a,  b  and  g ).  Rather  than  doing  this,  we  assume  that  solutions  of  (2.1 
-  2.5)  or  (2.8  -  2.11)  satisfy  certain  a-priori  estimates,  namely: 

For  fixed  T  >  0  there  are  positive  constants  n  and  M,  depending  on  norms  of  the  initial 
data  and  T ,  such  that  any  classical  solution  (i?(x,  t),  u(x,  £))  of  (2.1  -  2.5)  on  Qt  satisfies 

|u>(x,f)|  <  M  ,  |u(x,t)|  <  M  (2.15) 

and 

<pq(x,u(x,t),vx(x,t))>  n>  0,  (2.16) 

for  (x,  t)  e  Qt\  correspondingly,  if  (w(x  ,t),u(x,t))  is  a  classical  solution  of  (2.S  -  2.11)  on 
Qt,  then  (2.15)  and 

a(x,u(x,t),w(x,t))  >  fi  >  0  (2.17) 

hold  for  (x,  f)  e  Qt- 

The  objective  is  to  reveal  (2.15  -  2.17)  as  a  “minimal’1  set  of  a-priori  estimates  sufficient  for 
continuation  of  solutions  in  some  appropriate  function  classes.  Although  uniform  parabol- 
icity,  embodied  in  (2.16)  or  (2.17),  is  not  in  general  necessary  for  well-posedness,  in  light  of 
the  phenomena  under  consideration  and  for  technical  simplicity  solutions  will  be  continued 
up  to  the  first  time  that  uniform  parabolicity  fails.  For  the  models  at  hand  (2.15  -  2.17) 
are  established  in  Sections  3  and  4.  Finally,  it  is  shown  in  Lemma  2.3  that,  under  natural 
restrictions  on  the  initial  data,  (2.15  -  2.17)  always  hold  provided  T  is  sufficiently  small. 

The  first  goal  is  to  prove  an  existence  theorem  for  (2.8  -  2.11).  For  this  the  initial 
data  are  taken  smooth 

w0(x)  e  C2+Q[0, 1] ,  u0(z)  e  [CQ[0, 1]] N  ,  (2.18) 

for  some  0  <  a  <  1,  and  compatible  with  the  boundary  conditions: 

wor(0  =  0  ,  i  =  0, 1 ,  (2.19)v 

in  case  (2.10)v  applies,  or 

w0(i)  =  0,  a(i,uo(i),0)w0zx(i)  +  6(i,uo(i),0)  =  0 ,  i  =  0, 1 ,  (2.19)s 

in  case  (2.10)5  applies.  We  prove: 

Theorem  2.1.  Let  (iv0(x),u0(x))  satisfy  (2.18),  (2.19)  and  assume  that  the  a-priori 
estimates  (2.15)  and  (2.17)  hold  for  some  T  >  0,  with  M  and  n  positive  constants 
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depending  at  most  on  ||u>o||2+a,  ||uo||a  and  T.  There  exists  a  unique  solution 
( w(x,t),u(x,t ))  of  (2.S  -  2.11)  on  Qt  such  that  w,  wt,  wx,  wxx  are  in  C'°f>°72(Q7’) 

and  u ,  ut  are  in  [CQ’Q/2((5r)]  N. 

Proof.  The  proof  of  uniqueness  is  lengthy  but  routine  and  it  is  omitted. 

In  view  of  (2.15)  and  (2.17),  the  triplet  (x,  u(x,  t),  w(x,  t))  takes  values  in  the  set 
E  =  {(x,  u,  w )  e  [0, 1]  x  Rn  x  R  :  |u|  <  M,  |u>|  <  M,  a(x,  u  ,w)  >  p).  By  modifying, 
if  necessary,  the  functions  a,  b  and  g  outside  some  open  set  containing  F,  it  is  assumed 
for  the  existence  part  of  the  proof  that  all  the  functions  involved  are  bounded,  globally 
Lipshitz  and,  wherever  appropriate,  with  globally  Lipshitz  derivatives.  Moreover,  that 

a(x,p,r)  >  |  >  0  (2.20) 

for  (x,p,  r)  e  [0, 1]  x  RN  x  R.  All  bounds  and  Lipshitz  constants  depend  only  on  M  and  p. 
For  the  remainder  of  the  proof  K  will  stand  for  a  generic  constant  that  can  be  estimated 
solely  in  terms  of  M,  p  and  T. 

We  work  with  the  boundary  conditions  (2.10)s;  the  boundary  conditions  (2.10)^  are 
treated  similarly.  Let  B  denote  the  Banach  space 

G  =  {w(x,t)  eCd'^2(QT)  :  w(0,t)  =  u>(l,t)  =  0,  0  <  t  <  T}  (2.21) 

and  let  C  stand  for  the  closed  subset  of  [C^/2(Qt)]  N 

C  =  {u(x,  t)  e  [C0^/2(Qt)}  N  :  u(x,  0)  =  u0(x),  0  <  x  <  1}  .  (2.22) 

For  our  purposes  (3  =  min{a,\}.  Define  the  map  T  :  B  —>  C  that  carries  W(x,t)e£>  to 
U(x,t)  the  solution  of  the  family  of  initial  value  problems 

Ut  =  g(x,U,  W(x,t))  ,  0  <  X  <  1,  0<i<T,  03 

U(x,0)  =  uo(x )  ,  0<x<l. 

Also,  for  A  e  [0, 1]  define  a  second  map  S  :  [0, 1]  x  B  xC  — *  B  which  takes  (A,  W{x,  t),  U{x,t)) 
to  w(x,t )  the  solution  of  the  initial-boundary  value  problem 

wt  =  a\(x,  U(x,t),  IF(x,  t))wxx  +  A  6(x,  U(x,t),  W(x,t ))  +  (1  -  A)F(x) , 

w(0,t)  =  w(l,t)  =  0  ,  0  <  t  <  T ,  (2.24) 

ie(x,0)  =  w o(x)  ,  0  <  x  <  1 , 

where  a\(x,p,r)  :=  Aa(x,p,  r)  -f  (1  -  A)^  >  by  (2.20),  and  F(x)  :=  w0zz(x ). 

Given  the  maps  T  and  <5,  construct  the  composite  map  P  :  [0, 1]  x  B  — ►  B  which 
carries  A  e  [0, 1],  W  eB  to 


w  =  P(A,  IV)  :=  S{\,  W,  T(IV)) . 


(2.25) 
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Observe  that  if  w(x,t)eB  is  a  fixed  point  of  P(l,  •)  and  u(x,t)  the  corresponding  solution 
of  (2.23),  then  ( w(x ,  t ),  u(x,  t ))  satisfies  (2.8  -  2.11)  on  Qt ■  Our  objective  is  to  demonstrate 
that  the  map  P  fulfills  the  hypotheses  of  the  Leray-Schauder  fixed  point  theorem  (for  a 
formulation  see  [10],  also  [6]).  To  this  end,  certain  properties  of  the  maps  T  and  S  are 
recorded  below. 

First  consider  the  map  T  :  B  —*  C.  Since  g  is  bounded  and  globally  Lipshitz,  the 
standard  theory  of  ordinary  differential  equations  implies  that  given  any  W(x,t)  eB  there 
is  a  unique  solution  U(x,t )  of  (2.23)  defined  on  [0, 1]  x  [0,  T]  and  such  that 

\U(x,t)\  +  \Ut(x,t)\  <  Iu.  (2.26) 

Moreover,  by  first  integrating  (2.23)i  for  two  distinct  points  xi,  X2  in  [0, 1]  over  [0,#],  0  < 
t  <T,  and  then  estimating  the  difference  using  Gronwall’s  inequality,  we  deduce  with  the 
help  of  (2.26) 

WUWe.Mi  <  Kiduola  +  Wh.tft  +  1)  •  (2.27) 

Next,  consider  the  map  S  :  [0, 1 }  x  B  x  C  — ►  B.  The  classical  Schauder  theory  for 
parabolic  equations  [8,  9]  implies  that,  given  any  triplet  (A,  W,  U)  e  [0, 1]  x  B  x  C,  there  is 
a  unique  solution  w(x,t)  of  (2.24)  on  Qt  belonging  to  C2+f3,1+^/2(Qr)  and  satisfying 

IMta  +  0,0/2  ^  Ai(|wo|2+a  +  ||W||/M/2  +  \\U\\i)J/2  +  1)  •  (2.28) 

The  constant  Ai  can  be  estimated  solely  in  terms  of  ||Wj|/j,0/2>  ||Pj|0,0/2>  Af  and  T. 
Finally,  consider  the  map  P  :  [0, 1]  x  B  — *  B.  P  is  well  defined  by  (2.25).  Also: 

(i)  For  any  fixed  Ae[0,l],  P( A,-)  :  B  — ►  B  is  compact  and  continuous. 

Since  the  injection  C2+^'1+^2(Qt)  — ►  C^,i3^2(Qt)  is  compact,  (2.27)  and  (2.2S)  imply 
that  P(A,  •)  is  a  compact  map.  Let  {JFn}  be  a  convergent  sequence  in  B,  Wn  — +  W  in 
CP,i3/2(Qt)',  consider  wn  =  P( A,  Wn).  Since  P( A,  •)  is  a  compact  map,  along  a  subsequence 
ivn  — >  w  in  C P'®!2  (in  fact  in  C2+/J  ,1+^  /2,  for  any  /?'  <  /3).  One  easily  shows  that 
w  =  P(A,  W).  Since  P( A,  •)  is  single  valued,  wn  — ►  w  along  the  whole  sequence  and  P( A,  •) 
is  continuous. 

(ii)  For  any  bounded  subset  1C  of  B,  the  family  of  maps  P(-,W)  :  [0,1]  — »  B, 

IV  e/C,  is  uniformly  equicontinuous . 

Let  K  be  a  bounded  subset  of  B.  Fix  W e)C.  For  U  =  T(W)  and  A,  p  in  [0, 1],  let  w\,  wp  be 
the  respective  solutions  of  (2.24).  Note  that  w\  =  P(A,1F),  wp  =  P{p,W).  The  difference 
w \  —  wp  satisfies  the  parabolic  equation 

( w\  -  wp)t  =  a\(x,U(x,t),W{x,t))(iu\  -  wp)x z 

+  (A  -  p)[a(x,U{x,t),W(x,t))wPix  +  b(x,U{x,t),W(x,t))  -  (^wPii  +  F(. r))]  , 

(2.29) 

with  boundary  conditions  (2.24)2  initial  condition  ( w\  —  iC/,)(x,0)  =  0.  The  Schauder 
estimates  imply  that 

\\w\  -  Wph+0,1+13/2  <  A2|A  -  /?||||tep||2+^,l  +  /?/2  +  +  \\U\\l3,P/2  +  lu’o|2  +  a  +  l}  , 

(2.30) 
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with  A2  a  constant  as  in  (2.28).  Combining  (2.27),  (2.28)  and  (2.30)  we  arrive  at 

\\w\  ~  wp\\2+fi,i+fi/2  <  A3|A  -  p\  (2.31) 

with  A3  depending  only  on  fC,  the  initial  data,  p,  M  and  T.  Thus  P(-,  W),  W  eK  is  a 
uniformly  equicontinuous  family  of  maps. 

(iii)  P(0,  •)  has  precisely  one  fixed  point  in  B. 

For  A  =  0,  (2.24)i  becomes  the  heat  equation,  (2.23)  and  (2.24)  decouple,  and  (2.24)  has 
a  unique  solution  in  B. 

(iv)  Any  fixed  point  in  B  of  P( A,-),  0  <  A  <  1,  is  contained  in  some  bounded 
subset  K,  of  B. 

Let  uieSbea  fixed  point  of  P(A,  ■).  Set  u  =  T(w).  Then  (2.27)  implies 
u  e  jc^/^Qr)]  i  and  (2.26)  now  reads 

|u(z,f)|  +  |ut(x,f)|  <  K3  ,  (x,t)eQr.  (2.32) 

By  (2.27)  and  (2.28),  w  e  C2+/3'1+^2(QT)  and  satisfies 

wt  =  a\(x,  u ,  w)wXI  +  A  b(x,u ,  w)  +  (1  -  A )F(x)  (2.33) 

on  Qt  with  boundary  and  initial  conditions  as  in  (2.24).  Since  b  is  bounded,  the  maximum 
principle  yields 

t«’(^,OI  ^  SUP  |wo(x)| -f  A'4  +  (1  -  A)T  sup  |P(r)|,  ( x,t)eQr .  (2.34) 

0<r<l  0<i<1 


Next,  we  multiply  (2.33)  by  integrate  by  parts  over  [0,1]  x  [0,f],  and  use  Schwarz’s 
inequality  to  deduce 

nw2tdxdr  +  f  u^(x,  t)dx  <  K5  [l  +  /  u;or(x)<^x  +  (1  —  A)2  sup  |F(r)|2]  =:  C2  . 
Jo  Jo  0<r<l 

(2.35) 

Using  (2.35)  we  obtain 

| w(x,  t )  -  w(y,  f)|  <  C\x  -  y|1/2  .  (2.36) 


Also,  for  fixed  6  >  0  the  calculus  inequality 


rx  +  t>  rX+6  ft 

26\w(x,t)  -  u;(z,r)|  <  /  | w(x,t)  -  w(y,t)\dy  +  /  |  /  w3(y,s)ds\dy 

J x—b  J z— 5  Jr 


+  /  |u>(y,r)  -  w(x, 

Jx-5 


(2.37) 


r)\dy 


holds.  We  estimate  (2.37),  using  (2.36)  and  (2.35),  and  in  the  resulting  inequality  we  set 
S  =  |f  —  r  I1/2  to  arrive  at 

|u;(r,<)  —  ie(r,r)|  <  3C\t  —  r\ 1/4  .  (2.3S) 
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On  account  of  (2.34),  (2.36)  and  (2.38),  any  fixed  point  weB  is  contained  in  a  bounded 
set  of  C1/2,1/4(0r)?  and?  since  (3  —  min{a ,  ^},  also  in  a  bounded  set  of  B. 

The  map  P  fulfills  the  hypotheses  of  the  Leray-Schauder  fixed  point  theorem.  Thus, 
the  map  P(l,  •)  has  a  fixed  point  in  B.  If  w  eB  is  such  a  fixed  point  and  u  =  T(u>),  then 
( w(x ,  t),  u(x,  <))  is  a  classical  solution  of  (2.8  -  2.11)  on  [0, 1]  x  [0,  T}.  | 

We  collect  in  Lemma  2.2  certain  a-priori  estimates  for  solutions  of  (2.8  -  2.11)  that 
serve  as  a  starting  point  to  develop  an  existence  theory  for  the  system  (2.1  -  2.5).  Estimates 
(2.41)  capture  the  regularizing  effect  of  the  parabolic  equation  (2.8). 

Lemma  2.2.  Let  (w(x,t),u(x,t))  be  a  classical  solution  of  (2.8  -  2.11)  on  Qt  satisfying 
(2.15)  and  (2.17).  Then 


IHU  1  <  Ci(l+  /  ^oz(x)dx)  , 

’4  Jo 

(2.39) 

\Wh,l3/2  <C  ' L  +  f  WQz(x)dx  +  |«o|«)  , 

Jo 

(2.40) 

where  /3  =  min{a,  j}.  Moreover,  for  any  x,  ye  [0,1],  s,re[t,T]  with  t  >  0 

\iv(x,r)  -  w(y,  t )|  <  ^|x  -  y|1/2 
|u>(x,r)  -  io(x,s)|  <  ~^\T  ~  sl1/4  • 


(2.41) 


The  constants  C i,  C2  and  C3  above  depend  only  on  p,  M  and  T. 

Proof.  Estimate  (2.39)  is  a  direct  consequence  of  (2.34)  and  (2.35)  with  \  =  1,  together 
with  (2.36)  and  (2.3S);  (2.40)  follows  by  combining  (2.27)  with  (2.39). 

To  show  (2.41),  first  multiply  (2.8)  by  a(zwu  w)  and  write  the  resulting  identity  in  the 

form 


3‘f 

j  W0(z) 


i 

a(x,u,£) 


d(  +  w2x 


=  (wwx)z  + 


wb(x,u ,  w) 
a(x,  u,  w) 


£au(x,u,0 
a2(x ,  u,  {) 


d(, .  (2.42) 


In  view  of  (2J7),  we  may  assume  that  (2.20)  holds.  Integrating  (2.42)  over  [0, 1]  x  [O.t], 
0  <  t  <  T,  and  using  (2.9),  (2.10),  (2.15)  and  (2.20)  we  obtain 


w2zdxdr  <  Ki  . 


(2.43) 


Next,  mv’tiply  (2.8)  by  and  integrate  by  parts  over  [0, 1]  x  [0, <];  by  estimating  the 

resulting  identity  via  (2.15),  (2.20),  Schwarz’s  inequality  and  (2.43)  we  conclude  that 


tu>2(.t,  r)dxdr  4-  t 


x ,  t)dx  <  K 2  • 


(2.44) 
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A'i  and  AT2  depend  only  on  p,  M  and  T.  The  derivation  of  (2.41)  from  (2.44)  is  similar  to 
the  derivation  of  (2.36)  and  (2.3S)  from  (2.35);  in  (2.41)  C3  =  \fKf.  | 

Lemma  2.3  guarantees  that  the  a-priori  estimates  (2.15)  and  (2.17)  required  in  the 
hypotheses  of  Theorem  2.1  are  always  valid  for  T  sufficiently  small. 

Lemma  2.3.  Let  (w(x,t),u(x,t))  be  a  classical  solution  of  (2.8  -  2.11)  defined  on 
[0, 1]  x  [0,  T*),  for  some  T*  >  0. 

(a)  Suppose  that 

mi  :=  min  a(x,  uo(x),  iy0(x))  >  0 .  (2.45) 

0<£<1 

Then  given  any  positive  constants  p,  M  with  p  <  mi,  there  is  T  <  T*  depending  on  p ,  M 
and  the  L°°-norm  of  w0xx  such  that 

\w(x,t)  -  w0(x)\  <  M ,  (2.46) 

\u(x,t)  —  u0(x)\  <  M  (2.47) 

and  (2.17)  hold  for  (x,t)eQr- 

(b)  Suppose  that,  in  addition , 

m2  :=  min  a(x,  u0(x),  r)  >  0 ,  (2.48) 

0<r<l  v 

wo-  <r<w 04. 


where  w o_  =  inf  w0(x) ,  te0+  =  sup  Wq(x). 

°<*<l  0<z<l 

0  <  p  <  m2,  M  >  0,  M _  <  wQ _  <  ii>0+  <  M+, 
p ,  M ,  M-  and  A7+  such  that  (2.17),  (2.47)  and 


Then  given  any  p,  M,  A/_,  M+  with 
there  is  T  <  T*  depending  solely  on 


M-  <  w{x,t )  <  M+ 


(2.49) 


hold  for  (x,t)  e  Qt- 

Proof.  Let  W  =  w  —  wo,  U  =  u  —  uq.  Then  (W(x,  t),  U (x,  <))  satisfy  on  [0, 1]  x  [0,  T*)  the 
differential  equations 


Wt  -  A{x,  U,  W)WIX  =  A(x,  U,  W)w0lx(x)  +  B(x,  U,  W)  (2.50) 

Ut  =  G(x,U,W)  (2.51) 

with  boundary  conditions  (2.10)  and  initial  conditions  W(x,0)  =  0,  17(x,0)  =  0,  for 
0  <  x  <  1;  the  functions  A,  B  and  G  relate  to  a ,  b  and  g  through  formulas  of  the  format 

A(x,U,W)  =  a{x,uQ(x)  +  U,w0{x)  +  W).  (2.52) 


First  consider  part  (a).  Since  a ,  wo  and  Uo  are  continuous,  (2.45)  and  (2.52)  imply 
that  there  is  p  >  0  so  that  if  0  <  x  <  1,  |(/|  <  p  and  |W|  <  p  then  ^4(x,  U ,  W)  >  p  >  0. 
Set  k  =  min{p ,  A/}.  To  complete  the  proof  of  part  (a)  it  suffices  to  show:  There  is  T  <  T* 
such  that  for  (x,t)eQ 7-  the  triplet 

(x,  U(x,  t),  W(x,  0)  e  Ei  :=  {(x,  U,  W)  e  [0, 1]  x  Rs  x  R  :  \W\  <  k,  \U\  <  k } .  (2.53) 
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Then  (2.17),  (2.46)  and  (2.47)  follow  from  (2.52)  and  (2.53). 

Clearly  (2.53)  holds  on  Qr  for  some  r  sufficiently  small.  Moreover  so  long  as  (2.53) 
holds,  the  maximum  principle  for  the  parabolic  equation  (2.50)  gives  the  bound 

\W(x,  t)|  <  Fit  (2.54) 

where  Fx  =  sup  |A(x,  U,  VF)u>0rx(x)  +  B(x,U,W) |.  Also,  if  Gx  =  sup  \G(x,  U,  W)|,  then 
E  i  Ei 

(2.51)  together  with  Gronwall’s  inequality  yield 


\U(x,t)\  <  Git.  (2.55) 

Finally,  (2.54)  and  (2.55)  imply  that  (2.53)  holds  on  Qt,  for  any  0  <  T  <  T*  with 

T  <  Since  Fx  depends  on  the  L°°-norm  of  woxx,  the  resulting  T  will 

exhibit  the  same  dependence.  Under  the  stronger  hypothesis  (2.48)  this  dependence  can 
be  avoided. 

Consider  now  part  (b).  By  virtue  of  (2.48)  and  the  continuity  of  a  and  u0,  there  are 
p,  p-,  p+  with  p  >  0  and  p _  <  wo-  <  wo+  <  p+  such  that  0  <  x  <  1,  \U\  <  p  and 
P-  <  w  <  P+  imply  a(x,u0(x)  +  U,  w)  >  p  >  0.  It  now  suffices  to  show  that  there  is 
T  <  T*  such  that  for  (x,f)eQT  the  triplet 

(x,  U(x,  f),  w(x,  t))eE2  :=  [0, 1  ]x{UeRN  :  \U\  <  k}  x  [jfe_,  k+] ;  (2.56) 


here  k  =  min{p,M},  k-  =  max{p_,Af_},  k+  =  min{p+,  M+}  and  <  w0-  <  u>0+  < 
k+. 

Let  B±  =  sup  max{0,  ±6(x,  u0(x)  +  U,  to)}  and  consider  the  comparison  functions 

Ei 

W±(x,t)  =  ±B±t  —  u>o(a0  +  u>o±-  On  account  of  (2.50),  so  long  as  (2.56)  holds,  the 
functions  IV±  satisfy  the  differential  inequalities 


dtW-  -  A(x,  U,  W)d*W-  <  dtW  -  A(x,  U,  W)d\W  <  dtW+  -  A(x,  U,  W)d2xW+ 
W-(x,  0)  <  W(x,  0)  =  0  <  W+(x,  0) 


(2.57) 


and,  by  (2.19),  corresponding  inequalities  at  the  boundaries.  Using  comparison  principles 
for  parabolic  equations  we  obtain 


—  B-t  -f  wo-  <  io(x,  t )  <  B+t  +  u>o+  • 


(2.58) 


Moreover,  (2.51)  yields 


\U{x,i) |  <  G2t , 


(2.59) 


where  (?2  =  sup  |#(x,  w0(x)  i-U,  t^)|.  Finally,  (2.58)  and  (2.59)  imply  that  (2.56)  holds  on 

Ei 

k+-w0+  u»p_  -k-  k  ' 

1+ 


Qt,  for  any  0  <  T  <  T*  with  T  <  win  {  k+  ,  U’°~B  ■  ,  | 
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Theorem  2.1  in  conjunction  with  Lemmas  2.2  and  2.3  give  rise  to  the  following  local 
existence  and  continuation  theorem  for  the  initial-boundary  value  problem  (2.8  -  2.11)  in 
Schauder  spaces. 

Theorem  2.4.  Let  w0(x)  e  C1+a[0, 1],  u0(x)  e  [Ca[0, 1]]  N  satisfy  the  compatibility  con¬ 
ditions  u;Oi(0)  =  u>oi(l)  =  0  in  case  (2.10)v  applies,  or  tuo(O)  =  Wq(1)  =  0  in  case 
(2.10)s  applies,  and  suppose  that  (2.4S)  holds.  Then  there  exists  a  unique  classical  solu¬ 
tion  (w(x,t),  u(x,  t))  of  (2.8  -  2.11)  defined  on  a  maximal  interval  of  existence  [0, 1]  x  [0,  T*) 
such  that,  for  any  0  <  r  <  T  <  T*,  w  is  in  CQ''a/2(Qr),  u,  ut  are  in  [Ca,Ql//2(Qr)]  jV  and 
wt ,  wx,  wxx  are  in  C“'“/!([0,l]x[r,r]).  In  case  T *  <  oo,  as  t  |  T*  , 

limsup  sup  (|u(x,f)(  -f  |u;(:r,f)|)  =  oo  (2.60) 

t\T-  0<x<l 

and/or 

liminf  inf  a(x,u(x,t),w(x,t))  =  0.  (2.61) 

tTT*  o<x<i  v  '  v  ”  v  ’ 

Furthermore,  if  w0(x)  e  C2+a[0, 1]  and  the  compatibility  conditions  (2.19)  hold,  then  wt. 
wx,  wxz  are  in  Ca,a//2(Qr)  for  any  T  <  T* .  If,  in  addition ,  u0(x)e  [C1+Of[0, 1]]  N ,  then 
ux ,  uxt  are  in  \C0‘'0‘12(Qt)\  A  for  any  T  <  T* . 

Proof.  We  work  with  the  boundary  conditions  (2.10)5;  the  case  of  (2.10)y  is  treated 
similarly. 

Let  iuo(r)fC1+0[0,l],  u0(x)e  [C" [0, 1]]  ^  be  given,  satisfying  the  compatibility  con¬ 
ditions  u-’o(O)  =  u'u(l)  =  0  as  well  as  (2.4S).  We  proceed  to  establish  a  local  existence 
theorem  for  (2.8  -  2.11).  First,  construct  approximating  sequences  {u.*on}  and  {u0n}  such 
that  ic0„(x)  and  u0„(x )  are  C°°-functions  on  [0,1]  satisfying  (2.19)5,  and,  as  n  -+  cc, 

tr0n  ->  wo  in  Cl[0, 1]  ,  |u)0„|i+o  <  A'lu'oli+c  (2.62) 

Uon  — ♦  Uo  in  [£[0,1]]"  ,  |u0„U<  A'|uoU  (2.63) 

with  K  a  fixed  positive  constant  (for  details  of  such  a  construction  see  [19]). 

Consider  the  problem  (2.8  -  2.11)  with  initial  data  (w0n(x),u0n(x)).  Referring  to  part 
(b)  of  Lemma  2.3,  let 

mn  =  mini  a(x,  u0„(.r),  r)  :  0  <  x  <  1  ,  inf  w0n(x)  <  r  <  sup  u>o„(.t)  >  .  (2.64) 

l  °<*<1  0<x<l  J 

Since  (wo(x),uq(x))  satisfy  (2.4S)  and  Uq(x)  is  continuous,  liminf  mn  >  0.  By  throwing 

7  n— *oo 

away  a  finite  number  of  terms,  if  needed,  we  may  assume  that  m0  :  =  infmn  >  0.  Fix 

n 

p,  M,  A/_  and  A/+  such  that  A/_  <  inf  inf  won(x)  <  sup  sup  u’on(-i’)  <  A f+, 

n  n  0<x<l 

M  >  0,  0  <  p  <  m0.  Theorem  2.1  in  conjunction  with  Lemma  2.3  imply  that  for  each 
n  =  1,2,...  there  is  a  classical  solution  (wn(x,t),un(x,t))  of  (2.8  -  2.11)  defined  on 
[0,1]  x  [0,  T„],  with  smoothness  as  in  Theorem  2.1,  and  corresponding  to  the  initial  data 
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(u>0„(x),uon(2))-  Moreover,  T0  :=  inf  T„  >  0,  and  on  the  domain  QTo  =  [0, 1]  x  [0,To]  the 
functions  (wn(x,t),un(x,t))  satisfy  the  uniform  bounds 

\un(x,t)  -  u0n(*)|  <  M  ,  M_  <  wn(x,t)  <  A/+  (2.65) 

and 

a(x,  u„(x,  t),  wn(x,  t))>  fx  >  0 .  (2.66) 

Using  (2.62),  (2.63),  (2.65)  and  (2.66),  relations  (2.39)  and  (2.40)  in  Lemma  2.2  together 
with  (2.9)  imply  that  on  Qt0 


\\wn\\ 0,0/2  <  K  >  \\un\\0,0/2  +  \\dtUn\\ 0J/2  <  K ,  (2.67) 

where  (3  =  mm{a,  and  I\  is  a  constant  independent  of  n. 

Since  the  injection  C^'^^2(Qt0)  C^' ^2(Qt0)  is  compact  for  /?'  <  (3,  (2.67)  implies 

that  there  are  subsequences  {uv}  and  {un'}»  as  well  as  functions  w(x,t)  and  u(x,t),  with 

u>eC^/2(0ro)  and  u,  dtue  \C^'^^(Qt0)]N ,  such  that 


Wn‘  — ►  w 
un>  — »  u  ,  dtun> 


dtu 


in  C?’i3''2(Qt0), 
in  [C^^2(QTo)]N  . 


(2.68) 


Clearly  (w(x,  t),  u(x,  t))  satisfies  (2.9),  (2.10)s  and  (2.11).  Using  (2.68)  together  with 
results  on  families  of  solutions  of  parabolic  equations  (cf.  Friedman  [8,  Sec.  3.6]),  it 
follows  that  ( w(x,t),u(x,t ))  is  a  classical  solution  of  (2.8).  The  stated  regularity  of  this 
solution  is  an  outcome  of  the  interior  and  boundary  parabolic  estimates  [8,  Sec.  4.7]. 
Uniqueness  follows  from  a  lengthy  but  routine  argument  that  is  omitted. 

If  u;o(x)  €  C2+"[0, 1]  and  satisfies  (2.19)s  then  Theorem  2.1  implies  that  wt ,  wx  and 
wxx  are  in  Co,a^2(QT0)-  Suppose  that,  in  addition,  u o(x)  e  [C1+o,[0, 1]]  .  Now  u(x,t)  sat¬ 

isfies  (2.9)  with  w ,  u>x  eCa’°/2(QT0)-  Using  standard  theorems  on  continuous  dependence 
for  ordinary  differential  equations,  together  with  estimates  in  the  spirit  of  the  derivation 
of  (2.-7),  leads  to  ux,  uxt  e  [C“’"/2(0ro)]  ^ ■ 

Finally,  Theorem  2.1  implies  that  the  solution  (w(x,  t),  u(x,  t))  can  be  continued  on  a 
maximal  interval  of  existence  [0,1]  x  (0,T*),  such  that  either  T*  =  oo,  or  at  least  one  of 
(2.60)  or  (2.61)  occurs.  | 

Next,  we  turn  to  the  initial-boundary  value  problem  (2.1  -2.5).  We  assume  that  the 
initial  data  satisfy 

v0(x)  e  C2+“[0, 1] ,  u0(x)  e  [C1+Q[0, 1]] N  ,  (2.69) 


m0  =  mini  </?9(x,  Uo(x),  ip(x,  Uo(x),  r)  :  0  <  x  <  1,  inf  u’o(r)  <  r  <  sup  u>o(r)>>0, 

(  °<*<1  0<r< 1  ) 

(2.70) 

where  u>o(z)  is  given  by  (2.6),  and  the  compatibility  conditions: 


u0(0)  =  t>o(l)  =  0  ,  u’ox(0)  =  te0r(l)  =  0  ,  (2.71)v 
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in  case  (2.4)  y  applies,  or 


t^o(0)  =  u>0(l)  =  0, 


(2.71)s 


in  case  (2.4)5  applies.  We  prove: 

Theorem  2.5.  Under  the  hypotheses  (2.69  -  2.71),  there  exists  a  unique  classical  solution 
(u(x,  <),  u(x,  <))  of  (2.1  -  2.5)  defined  on  a  maximal  interval  of  existence  [0,  i]  x  [0,  T*) 
such  that,  for  any  0  <  T  <  T*,  v,  vx,  vt,  vxx  are  in  C“’“/2(Qt)  an<L  ui  Ut,  ux,  uxt  are  in 
If  T-  <  co,  then,  as  t  f  T*,  at  least  one  out  of  (2.60)  or 

limmf  o^nf  i  <pq(x, n(x,  t ),  vx{x ,  <))  =  0  (2.72) 

occurs.  Finally,  w(x,t)  =  <p(x,u(x,t),vx(x,t))  satisfies  (2.41),  where  C3  only  depends  on 
mo  and  the  sup-norms  of  wo  and  Uo. 

Proof.  For  concreteness,  we  treat  the  boundary  conditions  (2.4)s.  Let  (t;0(x),  u0(x)) 
satisfying  (2.69),  (2.70)  and  (2.71)5  be  given  and  define  wo(x )  by  (2.6).  Consider  the 
problem  (2.8  -  2.11)  with  a,  b  and  g  defined  by  (2.12).  Theorem  2.4  asserts  that  there  is  a 
unique  solution  ( w(x ,  t ),  u(x,  t ))  of  (2.8  -  2.11)  defined  on  [0, 1]  x  [0,  T*)  and  with  regularity 
as  stated  there. 

Our  objective  is  to  define  v(x,t)  by  (2.13)  subject  to  the  initial  data  u(x.0)  = 
Vo(x),  0  <  x  <  1.  Then  (v(x,  t),  u(x,  t))  is  a  solution  of  (2.1  -  2.5).  For  v(x,t)  to  be 
well  defined  it  is  at  least  required  that,  for  each  fixed  xe[0, 1],  wx(x,-)  is  integrable.  In 
view  of  Theorem  2.4,  to  complete  the  proof  it  suffices  to  show  that,  for  some  r  small, 
wx ,  vt,  vxz  e  Ca'al2(Qr),  ux ,  uxte  [Ca'al2{Qr)\  N  and  also  that  (3.41)  holds.  This  is  ac¬ 
complished  by  a  density  argument. 

Consider  approximating  sequences  {u>on}  and  {u0n},  consisting  of  C°°-functions  on 
[0,1],  with  {u>0n}  as  in  (2.62)  and  {u0n}  satisfying 

Won  W0  m  [C'lO.l]]"  ,  lu0n  1 14-or  —  7iiluo|i+o  1  (2.73) 

as  n  — ►  oc.  Set  v0n(x)  =  v0(0)  +  d>(y,  u0n(y),  w0n(y))dy  and  observe  that 

v0n  — *  u0  in  C2[0,1]  ,  I^Onh+Q  —  -^(K’oh+a  +  |uoll-fa)  •  (2.74) 


Let  (wn(x,  t),  un(x,  t ))  and  (vn(x,t),  un(x,t ))  be  the  corresponding  solutions  of  (2.8  -  2.11) 
and  (2.1  -  2.5),  respectively;  vn(x,t)  is  defined  by  solving  (2.13)  subject  to  en(x,0)  = 
von(x).  Note  that  by  uniqueness  for  (2.8  -  2.11),  (2.6S)  and  (2.13)  imply  that 


wn  — ►  w  ,  dxvn  — >  dxv  in  C(Qt) 

un-*u  in  [C(Qt)]N 


for  any  fixed  T  <T*. 

The  functions  vn  and  u>n  :=  dxun  satisfy  the  equations: 

dtvn  =  An(x,t)dx vn  +  Bn(x,t)  ■  u>n  +  Cn(x,t ) 


(2.75) 


(2.76) 


15 


dtun  =  Dn(x,  t)un  +  En(x,  t)dlvn  +  F„(x,  t ) ,  (2.77) 

where,  on  account  of  (2. 66),  (2.67)  and  (2.13)i,  the  components  of  An  —  Fn  are  uniformly 
bounded  in  Ca,a^2(QT),  and  for  n  large  An(x,t)  >  (j.  for  some  /i  >  0.  The  Schauder 
estimates  for  (2.76),  together  with  (2.74),  yield 

||wn||2+a,l  +  a/2  <  Ks[\vo b+a  +  No|l+a  +  \\<^n\\a,a/2  +  1]  •  (2.78) 

Proceeding  as  in  the  derivation  of  (2.27),  (2.77)  implies  that  on  Qr  with  0  <  r  <  T 

IMka/2  <  iCi[|uo|i+a  +  (r  +  r1-“)||ai2un|U)a/2  +  1] .  (2.79) 

The  constants  K3  and  K\  are  independent  of  n.  Combining  (2.78)  and  (2.79)  we  conclude 
that,  provided  r  -f  r1-a  <  2k\k^  '  es^ma^e 

lkn(|2+<»,l+a/2  +  ||wn||a,o/2  <  F'sduoh  +  a  +  [«o|l+a  +  1)  (2.80) 

is  valid  on  Qr.  Also,  by  virtue  of  Lemma  2.2,  u>n(x,t)  satisfies  (2.41). 

Relations  (2. SO)  and  (2.75)  give  wx ,  vxx  €  Ca,a^2(Qr)  and  uxe  [C“,o,/2(0r)]  N.  Also, 

vt  —  wxe  Ca’a/2{Qr)  and  uxt  =  fx  +  fp  •  ux  -f  /,  •  vxx  e  [C'a,Qf/2(Qr)] N  ■  Finally,  u>(x,f) 
satisfies  (2.41).  g 


3.  On  the  Competition  of  Strain  Softening  and  Strain  Rate  Dependence 


The  scope  of  this  Section  is  to  elucidate  the  competition  between  the  destabilizing 
influence  of  strain  softening  and  the  stabilizing  influence  of  strain  rate  sensitivity  and  their 
effect  on  the  response  of  shearing  motions.  Also,  to  provide  quantitative  criteria  that 
determine  which  one  prevails. 

We  use  as  a  test  case,  the  initial-boundary  value  problem  consisting  of  (1.1  -  1.3), 
namely 

Vt  =  (r(u)vj). 

0  <  x  <  1,  t>0  (3.1) 


ut  =  vx 


with  boundary  conditions  (1.9)s  and  initial  conditions  (1.5).  The  initial  data  are  taken 
smooth:  t>0(z)  e  C2+o,(0, 1],  ito(z)  e  C1+o[0, 1],  for  some  0  <  a  <  1;  compatible  with  the 
boundary  conditions:  ctq(0)  =  a 0(1)  =  1;  and  satisfying  the  sign  restrictions 


u0(x)  >  0,  <7o(x)  >0  ,  0  <  x  <  1 . 


(3.2) 


Henceforth,  we  will  refer  to  this  problem  as  ( V)s •  Recall  that  r(u)  is  a  smooth  function 
satisfying  (1.4)  and  n  is  a  positive  parameter. 

The  theory  developed  in  Section  2  implies  the  existence  of  a  unique  solution 
(u(x,t),u(x,<))  of  ( V)s ,  defined  on  a  maximal  interval  of  existence  [0,1]  x  [0,r*),  such 
that  v,  vt,  vz,  vxx,  u,  ut  and  uz  are  in  Ca,a/2(QT),  for  any  T  <  T*.  In  addition,  given 
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any  compact  subset  fC  of  (0,  oo)  x  R,  (a(x,t),u(x,t))  escapes  K  as  t  |  T"\  i.e.,  there  are 
sequences  {x„}  C  [0,1]  and  {<„}  with  tn  |  T*  such  that  (er(xn,  tn),  u(xn,  tn))  g  K.  The 
identification  of  ( V)s  with  (2.1  -  2.6)  is  done  by  setting  w  =  ip(x,p,q)  =  T{p)qn  —  1. 
(Although  is  not  C2  at  q  =  0,  it  is  C 2  when  restricted  to  compact  subsets  of  [0, 1]  x  R  x 
(0,oo);  this  remark  together  with  the  results  of  Section  2  provide  the  above  statements). 

Our  objectives  are:  (a)  to  characterize  the  class  of  functions  r(u)  and  parameters 
n  that  guarantee  global  solvability  for  ("P)s,  and  (b)  to  study  the  behavior  of  solutions 
( v(x,t),u(x,t ))  of  ('P)s- 

To  this  end,  it  is  expedient  to  use  a  different  formulation  of  (3.1).  Note  that  for 
(x,t)  €  [0, 1]  x  [0,T*),  a  solution  of  ( V)s  satisfies 


cr(x,t )  >  0,  ut(x,t)  =  vx(x,t)  >  0,  u(x,f)  >  uo(z)  >  0.  (3.3) 


A  simple  calculation,  using  (1.1  -  1.3),  shows  that  (<x(x,  t),  u(x,  t))  is  a  positive  solution  of 
the  reaction-diffusion  system 


(crn)t  =  r(u)-azx  + 


1  t'(u ) 
n  r(u)1  +  n 


r 


(3.4) 


ut 


X 

(Tn 


r(u)n 


(3.5) 


on  [0, 1]  x  [0,T*).  For  this  step  and  for  the  remainder  of  the  Section  we  assume  that 
(u(x, f), u(x, f))  enjoys  some  additional  smoothness,  namely,  vxt ,  vxzx  and  urr  are  in 
Co,“/2(0r))  f°r  any  T  <  T* .  Such  solutions  axe  generated  if  we  take  smoother  initial  data, 
a  hypothesis  that  can  later  be  relaxed  using  density  arguments  (cf.  Section  2).  Integrating 
(3.5),  yields 

$(u(x,f))  =  $(u0(x))  +  f  a^(x,  r )dr  ,  (3.6) 

Jo 

where 

$(“)  =  (3-7) 

An  important  ingredient  of  the  forthcoming  analysis  lies  in  estimating  cr(x,t)  by  means 
of  comparison  principles  (e.g.  [13,  Ch.  3,  Sec.  7])  for  the  parabolic  equation  (3.4);  in  turn, 
u(x,f)  is  estimated  using  (3.6).  We  state  the  comparison  principle  used  as  a  lemma  for 
future  reference. 


Lemma  3.1.  Suppose  that  for  any  T  <  T*,  <71(x,f)  and  cr2(x,t)  are,  both  in  C2,1((5t)  wri/i 
<7j (x,  f )  >  0  and  <72(x,f)  >  0,  A(x,t)  £  C(Qt )  with  A{x,t )  >  0  and  B(x,t)  €  C(Qt)-  If- 
for  any  T  <  T*, 


(<7f)(  -  A(x,t)crlxx  +  B(x,t)(af  )2  <  (<rf)t  -  A(x,t)a2zz  +  B{x,t)(a£  )2,  on  QT, 

<Ti(i,t)  <  cr2(f,<)  ,  i  =  0, 1,  0<<<T, 

<7i(x,0)  <  a2(x,0)  ,  0<x<l 

(3-8) 
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then 

<r\  (x,t)  <  <T2{x,t)  ,  0<x<l,0<t<T*.  (3.9) 

In  applications  of  the  lemma,  the  functions  A  and  B  axe  taken  respectively,  A  =  r(u)  »  >  0 
and  B  -  ~n~r TTT  >  0. 

n  r(u)1+  n 

Our  first  theorem  characterizes  global  solvability  of  ( V)s  in  terms  of  the  behavior  of 
$(u)  as  u  — ►  oo. 

Theorem  3.2.  Let  (v(x,t),u(x,t))  be  a  classical  solution  of  ('P)s  defined  on  a  maximal 
interval  of  existence  [0,1]  x  [0,T*).  Suppose  that  r(u )  satisfies  (1.4).  Then: 

(i)  T*  =  oo  if  and  only  if  $(oo)  =  oo. 

(ii)  If  T *  <  oo,  then 

lim  sup  u(x,t)  =  oo .  (3.10) 

t—T *  0<x<l 


Proof.  Let  (v(x,  t),  u(x,  /))  be  a  solution  of  ('P)s  on  [0, 1]  x  [0,T*),  with  T*  maximal;  let 
cr(x,t)  be  defined  by  (1.3).  We  estimate  the  solution  in  the  interval  of  existence  [0,1]  x 

[o.n. 

Under  hypothesis  (1.4),  any  positive,  concave  function  E(x)  satisfies 

— r(u)n£ri  -  >  °-  (3.H) 

n  r(u)1+n 

If,  in  addition, 

E(i)  >  a0(x)  ,  0  <  x  <  1  (3.12) 

then  applying  Lemma  3.1,  with  comparison  functions  <r(x,t)  and  £(x),  yields 

<r(x,  t)  <  S(x)  .  (3.13) 


First,  we  show  (ii).  Assume  that  T*  <  oo  and  at  the  same  time  u(x,f )  is  bounded 
from  above  on  [0, 1]  x  [0,T*),  i.e., 


«o-  <  u(x,t)  <U+  <  oo; 


(3.14) 


here  we  used  the  fact  that  u(x,t)  >  inf  Kn(i)  =:  uo_.  Let  s(t)  be  the  solution  of  the 

_  0<r<I 

initial  value  problem 


—  S”  +  Bq(s  n  )2  =0 
s(0)  =  <7q—  :=  inf  cr0(x) , 

0<i<1 


where  Bq  =  max - - — Integrating  (3.15),  yields 

«o-<«<t/+  n  r(0  +" 


s(t)  = 


oo- 


(l  +  &o-Bot) 


>  0 


(3.15) 


(3.16) 
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On  account  of  (3.15)  and  (3.4),  the  comparison  functions  s(t)  and  cr(x,t)  satisfy  the 
parabolic  differential  inequality  (3.8)  on  Qt,  for  any  T  <  T*.  Using  Lemma  3.1,  we 
deduce 

s(t)  <a(x,t).  (3.17) 

Estimates  (3.12),  (3.14)  and  (3.17)  imply  that  (cr(x,  t),  u(x,  t))  remains  in  some  compact 
subset  of  (0,oo)  x  R  as  t  f  T*.  But  then  the  solution  can  be  continued  past  T*,  which 
contradicts  the  assumption  that  T*  is  maximal  and  finite.  Since  ut>  0,  we  conclude  that 
if  T*  <  oo  and  maximal  then  (3.10)  holds. 

Next,  we  proceed  to  prove  (i).  Consider  $(u)  defined  for  u  €  (0,  oo)  by  (3.7).  $(u) 
is  increasing  and  invertible  with  an  inverse  function  <£-1(£)  defined  for  £  €  ($(0),  $(oo)) 
and  increasing.  In  case  $(oo)  =  oo,  combining  (3.6)  and  (3.13),  we  arrive  at 

u(x,t)  <  (<3?(u0(x))  -f  £"(x )t)  .  (3.18) 

Then  (3.14)  holds  for  any  Qt  with  T  >  0  and,  necessarily,  if  T*  is  maximal  then  T*  =  oo. 
By  contrast,  in  case  $(oo)  <  oo,  (3.6)  at  x  =  0  or  x  =  1  together  with  (1.9)5  leads  to 

$(u(i,<))  =  $(u0(i))  +  t  ,  i=0,l.  (3.19) 

In  turn,  (3.19)  implies 

u(i,t)  — >  00  as  t  —*  Ti,  (3.20) 

where  T,  —  $(00)  —  4>(u0(i))  <  00  for  i  =  0,1.  Therefore,  in  case  $(00)  <  00,  T*  < 
mm{T0, Tj }  <  00.  | 

According  to  Theorem  3.2,  the  criterion  for  global  solvability  of  ( V)s  is  the  divergence 
of  the  integral  r(£)*  d£  =:  4>(oo).  Therefore,  it  is  the  decay  rate  of  r(u)  as  u  — ►  00 
that  determines  global  existence  for  ('P)s-  The  cleiss  of  positive,  decreasing  constitutive 
functions  t(u)  can  be  decomposed  into  two  categories,  depending  on  whether  4>(oo)  is 
finite  or  infinite.  Roughly  speaking,  the  dividing  line  consists  of  functions  r(u)  that  decay 
to  zero  like  the  power  u~n. 

Next,  we  restrict  attention  to  functions  r(u)  such  that  $(00)  =  00  (and  thus  T*  =  00) 
and  consider  the  asymptotic  behavior  of  solutions  (u(x,  <),  u(x,  t))  of  (V)s  as  t  — ►  00.  In 
case  r(u)  =  r0  a  constant,  a(x,t)  is  a  positive  solution  of 

(<7n)<  =  TqCTix  (3.21) 

subject  to  the  boundary  conditions  (1.9)s;  thus  <j(x,t)  — ♦  1  uniformly  in  x  €  [0,1]  as 
t  —*  00.  The  question  is  whether  this  behavior  persists  for  positive  and  decreasing  functions 
r(u). 

Two  representative  classes  of  functions  r(u)  are  considered:  Class  (HI)  consists  of 
functions  that  decay  to  a  positive  constant  r( 00)  at  a  rate  dominated  by  a  power,  i.e.,  for 
some  c  >  0  and  o  >  0 

t(u)  >  t(o o)  >0  ,  0  <  —t'(u)  <  ,  u  >  0.  (HI) 
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Class  (H2)  consists  functions  that  decay  to  zero  like  a  power,  up  to  first  order  derivatives, 
i.e.,  for  some  c  >  0  and  m  >  0 


11  c 

<  r(u)  <  —  ,  0  < 

^  _ ,  m  —  V  /  —  . .  m  7 


t'(u)  m 


c  ir 


u' 


r(u) 


u  >  0 . 


u 


(H  2) 


Note  that,  to  guarantee  global  existence  for  functions  of  class  (H2),  we  assume  0  <  m  <  n 
so  that  $(oo)  =  oo.  Also,  note  that  under  hypothesis  (Hi), 


..  $(«)  .  ,-L 

lim  -  =  r(oo)n  , 

u— oo  U 


while,  under  hypothesis  (H2)  for  0  <  m  <  n 


and  for  m  =  n 


We  prove: 


$(u)  .  .  $(u) 

Inn  sup  — ; — <  oo  ,  lim  inf >  0 , 

U— OO  u  »  U— oo  ul  -» 


$(u)  .  -  $(u) 

lim  sup  — -  <  oo  ,  lim  inf  — -  >  0 


inu 


u-»oo  £nu 


(3.22) 


(3.23) 

(3.24) 


Theorem  3.3.  Suppose  the  function  r(u)  belongs  to  either  the  class  (HI)  or  the  class  (H2) 
with  0  <  ~  <  ^.  Let  (v(x,t),u(x,t))  be  a  classical  solution  of  ( V)s  on  [0,1]  x  [0,oo), 
corresponding  to  initial  data  (vQ(x),u0(x))  with  u0(x)  >  0,  <70(x)  >  0  for  0  <  x  <  1  and 
<r0(0)  =  <70(1)  =  0.  Then,  for  any  choice  of  the  initial  data  in  case  0  <  n  <  2  and  under 
restrictions  for  the  data  that  are  outlined  below  in  case  n  >  2,  the  following  hold:  As 
t  — >  oo, 

a(x,t)  =  l  +  0(t~l}),  (3.25) 

s-^ds)  (3.26) 

and 

v(x,t)  =  /  v0{y)dy+  /  f  ^  d(dy  +  0(t~d+'f)  (3.27) 

Jo  Jo  Jy  r(u(f,<))" 

uniformly  on  [0,1].  In  case  (HI)  holds  0  =  a  >  0  and  7  =  0,  while,  in  case  (H2)  holds 
0<  3  =  <  1  and  7  =  -2-. 

Proof.  Let  (u(i,  t),  it(x,  t))  be  a  classical  solution  of  (V)s  defined  on  Q ^  :=  [0, 1]  x  [0, 00). 
Then  (3.3)  holds  and  (cr(x,  <),  u(x,  t))  satisfies  (3.4),  (3.5)  on  Q oo-  We  proceed  to  obtain 
an  initial  a-priori  estimate,  independent  of  t,  for  a(x,t)  using  comparison  principles.  In 
the  sequel  I\  will  stand  for  a  generic  constant  that  can  be  estimated  in  terms  of  the  initial 
data  and  properties  of  the  function  r(u). 

Set  -7(11)  to  be  the  quotient  of  the  coefficient  of  the  reaction  term  over  the  coefficient 
of  the  diffusion  term  in  (3.4) 


7(u)  =  -(  - 


n 


r'(u)]  1 

r(u)  r(u)2/n 


(3.28) 
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If  we  can  find  positive  functions  5(x )  and  s(x)  defined  on  [0, 1]  and  satisfying  for  ( x,t )  6 
Q oo  the  differential  inequalities 


-Sxx+7(«(*,<))S2/n  >0 

S(x)  >  <70(x) 

(3.29) 

and 

-■Sn+7(w(x,f))s2/n  <  0 
s(x)  <  <To(x)  , 

(3.30) 

respectively,  then  Lemma  3.1 

[  implies 

s(x)  <  a(x,t )  <  S(x) 

(3.31) 

for  ( x,t )  €  Qoo- 

Next,  we  examine  the  possibility  of  finding  such  functions  5(x)  and  s(x).  Observe 
that  under  hypothesis  (HI), 


r  r 

0  <  7 (u)  <  —  <  I<2  ,  u  >  inf  u0(x )  (3.32) 

v  ua  ~  0<x<l 

while,  under  hypothesis  (H2)  for  0  <  ~  < 

0  <  7(u)  <  ■  <  AT4  ,  u  >  inf  u0(x) .  (3.33) 

U1  z  n  0<jt<1 

Under  either  of  (HI)  or  (H2),  any  concave,  positive  function  S(x)  satisfies  (3.29)i.  More¬ 
over,  for  any  choice  of  the  initial  function  <70(x)  >  0  there  is  a  concave  function  S(x )  such 
that  S(x)  >  <r0(x),  0  <  x  <  1.  Therefore,  the  right  hand  inequality  of  (3.31)  is  valid  for 
all  values  of  the  parameters  and  choices  of  the  initial  data. 

Turn  now  to  (3.30).  If  a  positive  function  s(x)  satisfies 

—Six  +  70S "  <0  ,  0  <  X  <  1  ,  (3.34) 

where  70  =  I\2  in  case  (HI)  holds  or  70  =  A'4  in  case  (H2)  holds  (0  <  ^  <  |),  then  s(x) 
also  satisfies  (3.30)  i.  Consider  two  cases  n  <  2  and  n  >  2. 

(i)  If  n  <  2,  the  parametric  family  s0(x)  =  |(x2  + 1)  satisfies  the  inequality  (3.34)  provided 
a  >  0  and  Q"_1  <  Given  any  (<70(x),  u0(x))  both  positive  on  [0,1],  we  fix  70  (which 
depends  on  t«o(x))  and  choose  a  sufficiently  small  so  that  sQ(x)  <  cro(x)  and  (3.34)  is 
fulfilled  for  x  6  [0, 1].  Then  5a(x)  satisfies  (3.30),  and  the  left  hand  side  of  (3.31)  is 
established. 

(ii)  If  n  >  2,  a  function  s(x)  satisfying  (3.34)  and  (3.30)2  can  only  be  found  for  restricted 
choices  of  the  data  ((70(x),  u0(x)).  For  instance,  given  cr0(x)  find  the  largest  a  >  0  such  that 
sQ(x)  =  f(x2  +  1)  <  cr0(x).  With  this  a  fixed,  s0(x)  satisfies  (3.34)  provided  70  <  a1-"  . 
In  view  of  the  choice  of  70  and  (3.32),  (3.33),  this  imposes  a  restriction  on  u0(x)  and/or 
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the  function  t{u).  If  these  restrictions  are  satisfied  then  the  left  hand  inequality  of  (3.31) 
holds. 

Henceforth,  we  restrict  attention  to  the  cases  when  (3.31)  holds.  Using  (3.31),  (3.6) 
yields 

it<$(u(x,t))-$(u0(x))<  Kbi  ,  t  >  0.  (3.35) 

As 

Combining  (3.3)  and  (3.35)  with  (3.22)  or  (3.23)  we  conclude  that,  under  hypothesis  (HI), 

+  1)  <  u(x,t)  <  I<6(t  4-  1)  ,  (3.36) 

A6 

while,  under  hypothesis  (H2)  with  0  <  ^  <  5, 


-7*(t  +  1) n~m  <  u(x,t)  <  K7(t  +  l)n-m  . 
A  7 


In  view  of  (1.1)  and  (1.9)s  we  have  the  following  identity: 


1  d_ 

2  dt 


v2dx  + 


(Ttvxtdx  =  0 . 


(3.37) 


(3.38) 


Using  (3.3S),  together  with  (3.4)  and  (1.1), 

1  d  f1  9  f1  x  1  _x  1  f1  /  t'(u)xct1+  n 

2 dii  v‘dx+n l  = L 


(3.39) 


which,  together  with  Schwarz’s  inequality,  implies 

1  d  f1  9  ,  n  [l  ,  s  1  ,_x  1  ,  ^  1  y1  /r'(u)\2  <^1+"  j_ 

/  v2dx  +  -  t(u)"<7  ”vxtdx<—  (—— )  -— rdx. 

2  dt  JQ  2  Jo  2n  y0  rvu)  r(u)" 

Also,  the  calculus  identity 

MM)  =  /  crr(y,t)dy  +  /  /  O'  xx  (MK, 

./O  Jo  Jy 


(3.40) 


(3.41) 


together  with  (1.1),  (1.9)s  and  Schwarz’s  inequality  lead  to 

Vt(x,t)<  (  vlt(Z,t)d£. 

Jo 


(3.42) 


Suppose  first  that  (H2)  holds  and  0  <  Then,  combining  (3.31),  (3.37)  and 

the  inequalities  in  (H2)  with  (3.40)  and  (3.42)  we  arrive  at  the  differential  inequality 


d_ 

dt 


L 


1  1 

v]dx  +  —  (t  +  1)‘ 
A  8 


n  .  3m-2n 

v2dx  <  K$(t  +  1)  n"m  . 


(3.43) 
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If  0  <  ^  <  |,  then  integrating  (3.43)  we  obtain 


/  V2t(x,t)d.X  <  Kl0(t  +  ; 

Jo 


(3.44) 


however,  if  ®  then  (3.43)  does  not  provide  any  decay.  Relation  (3.44),  in  conjunction 
with  the  Poincare  inequality 


(<r(*,f)-l)2<  f  aid,  t)d( 
Jo 


(3.45) 


and  (1.1),  yields  (3.25)  when  (H2)  holds. 

When  (Hi)  holds,  the  same  sequence  of  steps  using  (3.36)  and  (HI)  in  the  place  of 
(3.37)  and  (H2)  lead  to  the  differential  inequality 


i  Jo  V*dX  +  7^7  Jo  V*dX  ~  Kl2^  +  ^  ^  ’ 


which,  once  integrated,  yields 


[  v*(x,t)dx  <  I\n(t  +  1)' 
Jo 


(3.46) 


(3.47) 


Combining  (3.45)  and  (3.47)  we  arrive  at  (3.25),  in  case  (Hi)  holds. 

To  show  (3.26),  observe  that  (3.6),  (3.7),  Poincare’s  inequality,  (3.31)  and  (1.1)  yield 


/  T(0n^~  t  5:  /  |<7 "  (x, r)  —  l|dr 

Juo(z)  *  Jo 


<  A'14  f  (  f  Vf  (x,  r)dx)  1^2d.T  . 
Jo  Jo 


(3.4S) 


Using  (3.4S)  together  with  (3.44)  or  (3.47)  in  cases  (H2)  or  (HI),  respectively,  we  deduce 
(3.26). 

Finally,  the  identities 


v{x,t)=  f  v(y,t)dy  +  l  f  °  JMLd^dy 
Jo  Jo  jy  r(u(£,t))n 

and  i 

/  v(y,t)dy  =  /  vo(y)dy 
Jo  Jo 

(by  (1.1)  and  (1.9)s),  together  with  Poincare’s  inequality  and  (3.31)  imply 

„(*,<)-/  voMdy-l  £ 

vfdx)1/2(  f  —^dx). 

Jo  Jo  T(u)" 


(3.49) 


(3.50) 


(3.51) 
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Using  (3.47)  in  case  (HI)  holds,  or  (3.44)  and  (3.37)  in  case  (H2)  holds  we  arrive  at  (3.27). 

■ 

Between  the  class  of  functions  r(u)  satisfying  <£(oo)  <  oo,  for  which  solutions  of  ( V)s 
blow  up  in  finite  time,  and  the  class  of  r(u)  satisfying  (H2)  with  0  <  ^  for  which 
solutions  of  (V)s  behave  asymptotically  like  in  (3.21)  (cf.  (3.25)),  there  remains  a  gap  to 
be  analyzed  For  a  power  law 

t(u)  =  u~m  (3.52) 

the  gap  corresponds  to  the  powers  ^  <  1. 

In  the  sequel  we  consider  the  power  law  (3.52)  with  ^  and  analyze  the  behavior 
of  solutions  of  ('P)s-  Let  (v(x,t),u(x,t))  be  such  a  solution  defined  on  [0, 1]  x  [0,  T*);  here 
T*  =  +oo  if  \  <  <  1  and  T*  <  +oo  if  ^  >  1.  The  function  cr(x,t )  satisfies  (3.4)  with 

r(u)  as  in  (3.52),  while  (3.6)  yields 

it(x,i)1-^'  =  uo(aj)1-^’  +  (1  —  — -)  [  cr”(x,T)dr,  (3.53) 

n  J  o 

in  case  ^  ^  1,  and 

£nu(x,t)  =  £nu0(x)  +  (  <t*(x,t)cIt  (3.54) 

Jo 

in  case  —  =  1. 

We  examine  the  class  of  solutions  of  the  differential  inequality  (3.29).  Note  that,  by 
virtue  of  (3.28),  (3.52)  and  (3.3),  for  ~  \ 

7  (u(x,t))  =  —  u(x,t)2^-1  >  —  uo(x)2^"-1 .  (3.55) 

n  n 


Then,  Lemma  3.1  implies: 

Lemma  3.4.  Let  -^  >  j.  If  S(x )  is  a  smooth,  positive  function  satisfying  for  x  G  [0,1] 


-Sxx(x)  4-  7o(x)S"  (x)  >  0 
5(x)  >  cr0(x ) , 


(3.56) 


where  7o(x)  :=  -^uo(x)2^  *,  then 

a{x,t)  <  S(x)  ,  0  <  x  <  1,  0  <  t  <  T*  .  (3.57) 

If  7T  >  2»  glven  any  S(x)  >  <70(x)  >  0,  (3.56)  can  be  always  satisfied  by  choosing 
tto(x)  >  0  appropriately.  Namely  , 

m  ,  \ 2 m  —  l  Sxx(x)  ( 

-no(x)  -  >  TT7T-  (3-5S) 

n  S"(x) 

In  particular,  a 0(x)  can  be  used  as  a  choice  for  S(x)  for  restricted  choices  of  uQ(x). 
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(3.59) 


If  IT  =  2’  t;Iien  ^(x)  —  ^ofa)  >  0  will  satisfy  (3.56)  provided 

1  >  Sxr(x) 

2  ~  Si(x)  ' 

In  light  of  the  above  remark,  we  study  the  family  of  potential  choices 


SQ(x )  =  1  +  ai(i  —  1)  ,  0  <  x  <  1  (3.60) 

where  0  <  a  <  4.  Note  that  Sa(rc)  is  convex  and  attains  its  minimum  5a(|)  =  1  —  >  0. 

Also  5a(x)  — +  1  uniformly  on  [0, 1]  as  a  — ►  0. 

(a)  Case  ^  >  § 

If  the  initial  data  (uq( x),ao(x))  are  restricted  so  as  to  satisfy 


<70(x)  <  SQ(x)  ,  u0(x)  >  ua  := 


m 


(i-f)- 


,  0  <  x  < 


then  Lemma  3.4,  in  conjunction  with  (3.58),  implies 


a(x,t)  <  SQ(x)  ,  0<x<l,0<t<T*. 


(3.62) 


Observe  that  ua  — >  0  as  a  — »  0. 

Consider  now  the  problem  ( V)s  with  initial  data  (<70(a;),  Uq{x))  —  ( Sa(x),ua ),  for 
some  0  <  a  <  4.  Let  ( crQ(x,t),ua(x,t ))  be  the  corresponding  solution;  it  is  defined  on 
[0,1]  x  [0,T*).  Also, 

<r0(z,  t )  <  SQ(x)  =  1  +  ax(x  —  1)  (3.63) 

for  (x,t)  €  [0,1]  x  [0,T*). 

We  consider  three  separate  regions: 

(i)  2  <  n1  <  1:  Here  =  +°°.  In  addition,  (3.63)  and  (3.53)  yield 


ua(x,t)  < 


x_ia  , 

Ua  "  +  (1 - )Sq  ( x)t 

n 


(3.64) 


for  0  <  x  <  1,  0  <  t  <  oo.  Note  that,  since  5o(0)  =  5a(l)  =  1,  the  boundary  condition 
(1.9)5  together  with  (3.53)  imply  that  (3.64)  is  in  fact  an  equality  at  x  =  0  and  x  =  1. 
A  comparison  of  (3.25)  and  (3.26)  with  (3.63)  and  (3.64)  reveals  the  drastic  difference  in 
the  behavior  of  solutions  across  the  parameter  values  ^  In  particular  aQ(x,t)  does 

not  converge  to  1  as  t  — *  00,  and  spatial  nonuniformities  of  the  strain  u(x,t )  develop  and 
persist  in  time.  This  is  the  case  no  matter  how  close  to  the  constant  function  1  the  initial 
state  Sa(x)  is.  The  diffusion  is  in  this  case  too  weak  to  uniformize  the  solution. 

(ii)  ^  =  1:  The  situation  is  similar  to  part  (i)  with  (3.64)  replaced  by 

uQ(x,t)  <  uQexp^Sa  (z)tj  (3.65) 
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by  virtue  of  (3.54). 

(iii)  ®  >  1:  Here  T*  <  +oo.  The  solution  ( cra(x,t),ua(x,t ))  satisfies  on  [0, 1]  x  [0,T*)  the 
bounds  (3.63)  and 


ua(x,f)  < 


i-“ 

ua 


i  m  ,  x.  . 

~  |l  -  —  |S<J  (x)* 


(3.66) 


The  bound  on  the  right  hand  side  of  (3.66)  blows  up  for  the  first  time  at  the  boundary 
points  x  =  0  and  x  =  1  as  t  — ►  Tcr  :=  j — n  .  Then  (3.10)  together  with  (3.66)  imply 

that  T*  >  Tcr.  However,  since  Sa(0)  =  Sa(l)  =  1,  (3.66)  is  in  fact  an  equality  at  x  =  0 
and  x  =  1  and  thus  T*  =  Tcr.  The  function  uQ(x,t)  blows  up  exactly  at  the  boundary 
points  x  =  0  and  x  =  1  as  t  — *  Tcr  =  T*\  in  blowing  up  it  satisfies  (3.66)  and  appears  like 
two  shear  bands  located  at  the  boundaries  x  =  0  and  x  =  1.  Moreover,  as  t  — ►  T*,  the 
function  cra(x,<)  obeys  the  bound  (3.63),  while 


dtuQ(x,t)  =  dxva(x,t)  <  Ss  (x)  ul  "  —  |l - |Sa(x)t  m  "  (3.67) 

n  . 


with  an  equality  at  x  =  0  and  x  =  1. 

Note  that  in  the  above  cases  any  0  <  a  <  4  can  be  chosen.  Also,  by  choosing  other 
types  of  functions  S(x)  in  the  place  of  (3.60),  the  nonuniformities  that  develop  near  the 
boundary  can  be  made  very  strong,  at  the  expense  of  restrictions  on  the  initial  data  (cf. 
(3.58)). 

(b)  Case  “  =  5 

2. 

If  a  small  enough  so  that  2 a  <  |(l  —  7)",  then  ^(x)  satisfies  (3.59).  Consider 
initial  data  (SQ(x),  u0(x)),  with  a  small  and  uo(x)  >  0  but  otherwise  unrestricted.  Let 
(<70(x,  t),  ua(x,  t))  be  the  corresponding  solution  of  (7^)^;  it  is  defined  on  [0,1]  x  [0,oo). 
On  account  of  (3.60)  and  Lemma  3.4,  aQ(x,<)  satisfies  (3.63)  and  the  uniform  state  a  =  1 
again  loses  stability,  as  is  the  situation  in  the  case  ^ 


4.  The  Power  Law  and  Scale  Invariance 


For  the  particular  choice  t(u )  =  ~,  the  constitutive  relation  (1.3)  takes  the  form  of 
the  power  law 

<7  =  —  v?  (4.1) 

um 

where  m,  n  are  positive  parameters.  Under  (4.1),  (1.1  -  1.2)  read 

1 


t,(  =  (—v”) 

'  um  x  x 
ut  =  vt 

and,  correspondingly,  (3.4  -  3.5)  take  the  form 


(4.2) 


/  JL\  -UL  m  EL- 1  i 

(cr")  =u  "  azz - U"  o  0 

v  '  ‘  n 


m.  X 

Ut  =  U  n  <7  n  . 


(4.3) 
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The  power  structure  of  the  system  (4.2)  induces  the  following  scaling  property:  If 
(u(x,  t),  u(x,  t ))  is  a  solution  of  (4.2)  on  R  x  (0,  oo)  then  (ua(z,  t),  u\(x,  <)),  defined  by 

v\(x,  t)  =  A  “  u(Ax,  A ~°t)  (4.4) 

a  -^6  ■4*  1  | 

ua(x,<)  =  A  “  u(Ax,  A-«<),  (4.5) 

where  A  >  0  and  8,  a  are  any  constants  with 

—m(a  +  <5  +  1)  +  n(a  +  8)  +  a  —  8  +  1  =  0,  (4.6) 

is  also  a  solution  of  (4.2)  on  R  x  (0,  oo)  ([16]).  As  a  consequence 

a\(x,t)  =  A  s  er(Ax,A"«t),  (4.7) 

and  (cr\(x,t),u\(x,t))  satisfies  (4.3)  on  R  x  (0,oo)  for  any  A  >  0.  Therefore  the  systems 
(4.2)  and  (4.3)  are  invariant  under  the  group  of  stretching  transformations  T0is: 


Ax ,  t 


A -it, 


,  a 

A“  v 


a+S  +  l 


u 


u ,  a 


A  “  <7 


0  <  A  <  co 


(4.8) 


with  o,  8  constrained  by  (4.6). 

The  system  (4.2)  admits  a  special  class  of  solutions  describing  uniform  shearing 


u(x,f)  =  x  ,  u(x,t)  =  t  +  u0.  (4-9) 

They  correspond  to  initial  data  uo(x)  =  x  and  Cto(x)  =  uo,  where  u0  is  an  arbitrary 
positive  constant.  For  the  special  choice  Uq  =  0,  (x,f)  is  a  self  similar  solution  under  the 
transformation  TQ)$  with  a  =  —8  =  m~* . 

Consider  now  the  initial-boundary  value  problem  consisting  of  (4.2)  on  [0, 1]  x  {<  >  0} 
with  boundary  conditions  (1.9)y  and  initial  conditions  (1.5).  Suppose  that  the  initial  data 
are  smooth  v0(x )  6  C2+Q[0, 1],  uo(x)  €  C1+a[0, 1],  for  some  0  <  a  <  1,  they  axe  compatible 
with  the  boundary  data,  and  satisfy  the  sign  restrictions  <7o(x)  >  0,  u0(x)  >  0,  0  <  x  <  1. 
We  will  refer  to  this  problem  as  (V)v  (including  the  assumptions  on  the  initial  data). 

The  existence  theory  developed  in  Section  2  implies  that  {V)v  admits  a  unique  clas¬ 
sical  solution  defined  on  a  maximal  interval  of  existence  [0,1]  x  [0,T*).  Moreover,  if 
T*  <  +oo,  given  any  compact  subset  K  of  (0,  oo)  x  (0,oo),  (cr(x,t),u(x,t))  escapes  fC  as 
t  t  T*.  Also,  for  (x,t)  6  [0,1]  x  [0,T*), 

a(x,t)>  0  ,  ut(x,t)  >  0  ,  u(x,t)  >  u0(x) .  (4.10) 

The  uniform  shearing  solution  (4.9)  is  a  special  solution  of  (V)v  for  initial  data  (x,  uo). 
Our  objective  is  to  study  the  stability  of  this  solution.  To  this  end  we  use  a  transformation, 
motivated  by  the  scaling  properties  of  (4.2),  to  obtain  a  system  that  admits  invariant  re¬ 
gions  [4].  A  similar  idea  has  been  independently  pursued  by  Bertsch,  Peletier  and  Verduyn 
Lunel  [1]  for  a  related  system. 

We  prove 
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Theorem  4.1.  Suppose  m  <  n.  There  exists  a  unique  classical  solution  (v(x,  t),  u(x,  t)) 
of  (V)y  defined  on  [0,1]  x  [0,  oo)  such  that  v,  vz,  vt,  vxx,  u,  ux  and  ut  are  in  C0,'c‘/2(Qt), 
for  any  T  >  0.  Moreover,  if  m  <  min{n,  1},  as  t  — ►  oo, 


vx(x,t)  =  14-  0(t 13  x) 
u(x ,  t)  =  t  +  0(t P) 


and 


a(x,t)  =  t-m{l  +  0(tt-1)), 
uniformly  on  [0,1],  with  ft  =  mai{^,m}  <  1. 

Proof.  Introduce  the  transformations 

v(x,t)  =  V(x,s(t)) 

u(x,t)  =  (t  4-  1  )U(x,s(t)) 
a(x,  t)  =  (i  +  l)-mS(x,  s(f )) 


(4.11) 

(4.12) 

(4.13) 


with 


s(t)  =  tn(t  +  1)  , 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


which  are  motivated  by  the  form  of  the  uniform  shearing  solutions  and  the  scaling  invari¬ 
ance  (4.8).  Relations  (4.14  -  4.16)  and  (4.1)  induce 

£  =  ~;V: .  (418) 

Moreover,  since  (cr(x,  t),  u(x,  <))  satisfies  (4.3)  and 

<7*(0,0  =  <7x(1,<)  =  0  ?  oo,  (4.19) 

it  follows  that  (£(x,  s),  U(x,  s))  solve  the  system  of  reaction-diffusion  equations 

S.  =  n -  U'~*)  ^ 

V.  =  U*(Ei 

subject  to  boundary  conditions 

£r(0,s)  =  E*(l,s)  =  0  ,  s  >  0  (4.21) 


and  initial  conditions 

E(x,0)  =  cr0(x)  ,  f/(x,0)  =  u0(x)  ,  0  <  x  <  1 .  (4.22) 
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If  0  <  ^  <  1,  the  theory  of  Chueh,  Conley  and  Smoller  [4]  guarantees  that  (4.20) 
admits  positively  invariant  rectangles  of  arbitrary  size  in  the  first  quadrant  {(E,  U)  G  R2  : 
E  >  0,  U  >0}.  They  are  centered  around  the  line  E  =  Un~m  and  look  like  in  Fig.  1. 


E  =  Un~m 


Fig.  1:  Invariant  regions  for  (4.20) 


Given  initial  data  cr0(x )  >  0,  u0(x)  >  0,  let  U-,  U+,  E_  and  E+  be  the  defining  coordinates 
of  the  smallest  invariant  rectangle  containing  (er0(x),  uo(x)),  0  <  x  <  1.  Then 

E_  <  E(x,  s)  <  E+  ,  U-  <  U(x,s)  <  U+  .  (4.23) 

In  turn,  (4.15),  (4.16),  (4.14)  and  (4.18)  in  conjunction  with  (4.23)  yield 

U-(t  +  l)  <u(x,t)  <  U+(t  +  l),  (4.24) 

E_(t  4-  l)"m  <  <r(x,0  <  E+(i  +  l)"m  (4.25) 

E ZUT  <  vr(x,t)  <  .  (4.26) 

The  first  implication  of  (4.24  -  4.26)  is:  For  0  <  ^  <  1  the  functions  (<r(x,  <),  u(x,  t)) 
remain  in  a  compact  subset  of  (0,oo)  x  (0,oo)  for  any  finite  time.  Thus  T*  =  -foo, 
that  is  solutions  (v(x,  t),  u(x,  t))  of  ( V)v  are  globally  defined.  In  addition  (4.24  -  4.26) 
provide  preliminary  information  on  the  time  evolution  of  solutions.  They  are  supplemented 
below  with  parabolic-type  energy  estimates  to  establish  the  stated  asymptotic  behavior. 
In  what  follows  K  will  stand  for  a  generic  constant  that  depends  only  on  the  data  and  the 
parameters  m  and  n. 

Our  first  goal  is  to  estimate  the  £2-norm  of  Vf.  To  this  end,  differentiate  (1.1)  with 
respect  to  t  and  use  (4.1)  and  (1.2)  to  obtain 


vtt 


(-5- 

\um 


nv*  1 vxt  - 


m 

um+l 


,n+l 


(4.27) 
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We  multiply  (4.27)  by  vt  and  integrate  by  parts  over  [0,  1],  using  (1.9)y,  to  arrive  at 


v  —  m— 2 


Id  f1  f1  v?-1  2  ,  r1  u”+1  , 

2  dtj0  *  J0  u™  xt  J0  u "+1 

On  account  of  (4.24),  (4.26)  and  Schwarz’s  inequality,  (4.28)  yields 

Jt  J 1  vfdx  +  i(«  +  l)-"  jf’  vl ,dx  <  K,(t  +  1)' 

Finally,  combining  (4.29)  with  the  Poincare  inequality 

v2(x,t)<  j  v2xt(x,t)dx 
Jo 

we  arrive  at  the  differential  inequality 

[  v2dx  +  ~{t  +  l)~m  f  v2dx<Ki(t  +  iy 
dt  Jo  -hi  Jo 

Integrating  (4.31),  we  deduce 

[  v2(x,t)dx  <  (  f  v2(x,0)dx)exp{-jr-  /  (s  +  l)-mds) 

Jo  Jo  -*V  1  Jo 

+I< !  J  {s  +  l)~m-2exp{--^- J'{t +  \)~mdT}ds. 

In  case  m  <  1,  L’Hopital’s  rule  implies 

,.  fo(s  +  l)-m-2eip{-^-  f09(r  +  1  )~mdT}ds  _ 
hrn  - ; — - - -  =  A 1  . 

t-oo  (<  +  l)-2exp{7^-  /0(T  +  l)"mdT} 


v— m— 2 


In  view  of  (4.33),  (4.31)  yields  for  0  <  m  <  1 


V2(x,  t)dx  <  I\2(t  +  1)  2. 

By  contrast,  if  m  >  1,  (4.32)  does  not  provide  decay  for  the  £2-norm  of  vt- 
m  =  1  the  decay  rate  depends  on  the  coefficient  K\  in  (4.31). 

Equations  (4.3)2  and  (1.1)  readily  imply 


m  1  ft 

u'^(x,t)uz(x,t)  =  Uo^(x)uQx(x)+  -  /  a±-\x,T)vt(x,T)dT , 

n  Jo 


(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 
Finally,  if 

(4.35) 
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Two  cases  are  considered:  (i)  If  n  ^  1,  then,  by  virtue  of  (4.24),  (4.25)  and  (4.34),  (4.35) 
yields 


f1  \ux(x,t)\dx  <  K3(t  +  1)?  +  I<4(t  +  1)*  [\t  +  l)"“+m(  f1  v?(x,T)dx)1/2dr 
Jo  Jo  Jo 

<Ks(t  +  l)?  +IU(t  +  l)m. 

(ii)  If  n  =  1,  then  (4.35)  reads 

u~m(x,t)ux(x,t)  =  Uom(x)u0x(x)  +  v(x,t)  -  v0(x ).  (4.37) 


(4.36) 


Using  (4.24)  together  with  the  maximum  principle  for  (4.2)i  we  arrive  again  at  (4.36). 
Combining  (4.36)  with  the  identities 


u(x,t)~  /  u(y,t)dy  =  /  /  ux(Z,t)d£dy 

Jo  Jo  Jy 

/  u(y,i)dy  =  t+  U0(y)dy 
Jo  Jo 


(4.38) 


we  obtain  (4.12). 

Next,  use  the  identity 


nu"  1vxx  =  umvt  +  —i  ixv” 

u 

in  conjunction  with  (4.26),  (4.24),  (4.34)  and  (4.36)  to  deduce 

J  \vzx(x,t)\dx  <  I<7(t  +  l)m{J  v}(x,t)dx)X/2  +  J  |uI(x,t)|c?x 
<  A'^  +  ir^+A'^  +  l)^-1. 


(4.39) 


(4.40) 


Then  (4.11)  follows  by  virtue  of  (4.40)  and  the  Poincare  inequality 

K(x,*)-l|<  [  |urx(£,  t)\dx  .  (4-41) 

Jo 

Finally,  to  show  (4.13),  note  that  on  account  of  (4.11),  (4.12),  as  t  — >  oo 

<(r,t)  =  l  +  0(t/J"1)  (4.42) 

u-m(x,f)  =  rm(H-0(^"1)),  (4.43) 

where  =  mai(^,m}  <  1.  Combining  (4.1)  with  (4.42)  and  (4.43)  gives  (4.13).  | 

To  shed  some  light  on  the  relevance  of  the  constraint  m  <  1,  consider  the  case  1  < 
m  <  n  and  observe  that  (4.24)  implies 

U+m(t  +  l)“m  <  u~m(x,t)  <  UZm(t  +  l)"m  .  (4.44) 
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Thus,  the  diffusion  coefficient  in  (4.2)i  decays  like  (t  +  1)  m. 
Consider  the  problem 


vt  =  a(t)(v")x 


(4.45) 


on  [0,1]  x  (0,oo),  subject  to  (1.9)v  and  r(x,0)  =  uo(^)  with.  u0z(x)  >  0-  The  change  of 
variables 


V(x,s(t))  =  v(x,t) 

s(t)  =  f  a(r)dr 
Jo 

suggests  that  V(x,s)  satisfies 

V.  =  (■ vxn ), 


(4.46) 

(4.47) 

(4.48) 


subject  to  the  same  initial  and  boundary  conditions.  For  a(t)  =  (t  +  1)  m,  we  have: 
s0 o  :=  lim  s(t)  is  infinite  for  m  <  1,  but  finite  for  m  >  1.  Also, 

t— *oo 

limu(x,<)=  lim  y(x,s).  (4.49) 

t— «oo  a— -Joe 

If  Soo  =  +oo,  then  lim  u(x,<)  =  x;  however,  if  s oo  <  +oo,  in  general,  this  will  no  longer 
1—00 

be  true. 


References 

[1.]  M.  Bertsch,  L.  A.  Peletier  &  S.  M.  Verduyn  Lunel,  The  effect  of  temperature  dependent 
viscosity  on  shear  flow  of  incompressible  fluids,  SIAM  J.  Math.  Anal,  (to  appear). 

[2.]  N.  Charalambakis,  Stability  asymptotique  lorsque  t  —*  oo  en  thermoviscoelasticite  et 
thermoplasticite,  in  “Nonlinear  partial  differential  equations  and  their  applications: 
College  de  France  Seminar  Vol  IX”,  (H.  Brezis  and  J.-L.  Lions,  eds.),  Pitman  Res. 
Notes  in  Math.  No  181,  Longman  Scientific  &  Technical,  Essex,  England,  1988. 

[3.]  N.  Charalambakis  &  F.  Murat,  Weak  solutions  to  the  initial-boundary  value  problem 
for  the  shearing  of  non-homogeneous  thermoviscoplastic  materials.  Proc.  Royal  Soc. 
Edinburgh  Sect.  A  113  (1989),  257  -  265. 

[4.]  K.  Chueh,  C.  Conley  &  J.  Smoller,  Positively  invariant  regions  for  systems  of  nonlinear 
diffusion  equations,  Indiana  Univ.  Math.  J.  26  (1977),  372  -  411. 

[5.]  B.  D.  Coleman  &  M.  L.  Hodgdon,  On  shear  bands  in  ductile  materials.  Arch.  Rational 
Mech.  Anal.  90  (1985),  219  -247. 

[6.]  C.  M.  Dafermos  L.  Hsiao,  Global  smooth  thermomechanical  processes  in  one¬ 
dimensional  nonlinear  thermoviscoelasticity,  J.  Nonlinear  Analysis  6  (1982),  435  - 
454. 

[7.]  C.  M.  Dafermos  k  L.  Hsiao,  Adiabatic  shearing  of  incompressible  fluids  with  temper¬ 
ature  dependent  viscosity,  Quart.  Appl.  Math.  41  (1983),  45  -  58. 


32 


[8.]  A.  Friedman,  Partial  Differential  Equations  of  Parabolic  Type,  Prentice  Hall,  Engle¬ 
wood  Cliffs,  New  Jersey,  1964. 

[9.]  0.  A.  Ladyzenskaja,  Y.  A.  Solonnikov  &  N.  N.  Ural’ceva,  Linear  and  Quasilinear 
Equations  of  Parabolic  Type ,  (Translated  from  the  Russian  by  S.  Smith),  American 
Math.  Society,  Providence,  Rhode  Island,  1968. 

[10.]  J.  Leray  &  J.  Schauder,  Topologie  et  equation  fonctionelles,  Ann.  Sci.  Ecole  Norm. 
Sup.  51  (1934),  45  -  78. 

[11.]  A.  Needleman,  Material  rate  dependence  and  mesh  sensitivity  in  localization  problems, 
Comp.  Methods  in  Appl.  Mech.  and  Engineering  67  (1988),  69  -  85. 

[12.]  J.  A.  Nohel,  R.  L.  Pego  &  A.  E.  Tzavaras,  Stability  of  discontinuous  steady  states  in 
shearing  motions  of  a  non-Newtonian  fluid,  Proc.  Royal  Soc.  Edinburgh  Sect.  A  (to 
appear). 

[13.]  M.  H.  Protter  &  H.  F.  Weinberger,  Maximum  Principles  in  Differential  Equations, 
Prentice  Hall,  Englewood  Cliffs,  New  Jersey,  1967. 

[14.]  M.  Renardy,  W.  Hrusa  J.  Nohel,  Mathematical  Problems  in  Viscoelasticity,,  Pit¬ 
man  Monographs  and  Surveys  in  Pure  and  Applied  Mathematics,  Vol.  35,  Longman 
Scientific  &  Technical,  Essex,  England,  1987. 

[15.]  T.  G.  Shawki  &  R.  J.  Clifton,  Shear  band  formation  in  thermal  viscoplastic  materials. 
Mech.  Mater.  8  (1989),  13  -  43. 

[16.]  A.  E.  Tzavaras,  Ph.D.  Thesis,  Brown  Univ.,  Providence,  1985. 

[17.]  A.  E.  Tzavaras,  Shearing  of  materials  exhibiting  thermal  softening  or  temperature 
dependent  viscosity,  Quart.  Appl.  Math.  44  (1986),  1  -  12. 

[IS.]  A.  E.  Tzavaras,  Plastic  shearing  of  materials  exhibiting  strain  hardening  or  strain 
softening,  Arch.  Rational  Mech.  Analysis  94  (1986),  39  -  58. 

[19.]  A.  E.  Tzavaras,  Effect  of  thermal  softening  in  shearing  of  strain-rate  dependent  ma¬ 
terials,  Arch.  Rational  Mech.  Analysis  99  (1987),  349  -  374. 

[20.]  T.  W.  Wright  &:  J.  W.  Walter,  On  stress  collapse  in  adiabatic  shear  bands,  J.  Mech. 
Phys.  Solids  35  (19S8),  701-720. 

[21.]  F.  H.  Wu  &  L.  B.  Freund,  Deformation  trapping  due  to  thermoplastic  instability  in 
one-dimensional  wave  propagation,  J.  Mech.  Phys.  Solids  32  (1984),  119  -  132. 


33 


